Putting competing orders in their place near the Mott transition 
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We describe the localization transition of superfluids on two-dimensional lattices into commen- 
surate Mott insulators with average particle density p/q {p, q relatively prime integers) per lattice 
site. For bosons on the square lattice, we argue that the superfluid has at least q degenerate species 
of vortices which transform under a projective representation of the square lattice space group (a 
PSG). The formation of a single vortex condensate produces the Mott insulator, which is required 
by the PSG to have density wave order at wavelengths of q/n lattice sites (n integer) along the 
principle axes; such a second-order transition is forbidden in the Landau-Ginzburg- Wilson frame- 
work. We also discuss the superfluid-insulator transition in the direct boson representation, and find 
that an interpretation of the quantum criticality in terms of deconfined fractionalized bosons is only 
permitted at special values of q for which a permutative representation of the PSG exists. We argue 
(and demonstrate in detail in a companion paper: L. Balents et ai, cond-mat/0409470) that our 
results apply essentially unchanged to electronic systems with short-range pairing, with the PSG 
determined by the particle density of Cooper pairs. We also describe the effect of static impurities 
in the superfiuid: the impurities locally break the degeneracy between the q vortex species, and this 
induces density wave order near each vortex. We suggest that such a theory offers an appeafing 
rationale for the local density of states modulations observed by Hoffman et al, Science 295, 466 
(2002), in scanning tunnelling microscopy (STM) studies of the vortex lattice of Bi2Sr2CaCu20g+5, 
and allows a unified description of the nucleation of density wave order in zero and finite magnetic 
fields. We note signatures of our theory that may be tested by future STM experiments. 



I. INTRODUCTION 

One of the central debates in the study of cuprate su- 
perconductivity is on the nature of the electronic corre- 
lations in the 'underdoped' materials. At these low hole 
densities, the superconductivity is weak and appears only 
below a low critical temperature (Tc), if at all. Never- 
theless, there is no evidence of metallic behavior and an 
underlying Fermi surface of long-lived quasiparticles, as 
would be expected in a conventional BCS theory. This 
absence of Fermi liquid physics, and the many unex- 
plained experimental observations, constitute key open 
puzzles of the field. 

In very general and most basic terms, the passage from 
optimal doping to the underdoped region may be char- 
acterized by an evolution from conducting (even super- 
conducting) to insulating behavior. Indeed, most think- 
ing on the cuprates is informed by their proximity to a 
Mott insulating state, in which interaction-induced car- 
rier localization is the primary driving influence. Many 
recent attempts to describe the underdoped region have, 
by contrast, focused on the presence or absence of conven- 
tional orders (magnetism, stripe, superconductivity, etc.) 
as a means of characterization. This focus has practical 
merit, such orders being directly measured by existing ex- 
perimental probes. Theoretically, however, the immense 
panoply of competing potential orders reduces the pre- 
dictive power of this approach, which moreover obscures 
the basic conducting to insulating transition in action 
in these materials. In this paper, we describe a theo- 
retical approach in which Mott localization is "back in 



the driver's seat" . Less colloquially, we characterize the 
behavior near to quantum critical points (QCPs) whose 
dual "order parameter" primarily describes the emer- 
gence of an insulating state from a superconducting one. 
Remarkably, a proper quantum mechanical treatment of 
this Mott transition leads naturally to the appearance 
of competing conventional orders. The manner in which 
this occurs is described in detail below. 

Among recent experiments, two classes are especially 
noteworthy for the issues to be discussed in our paper. 
Measurements of thermoelectric transport in a number of 
cuprates^ have uncovered a large Nernst response for a 
significant range of temperatures above Tc. This response 
is much larger than would be expected in a Fermi liquid 
description. However, a model of vortex fluctuations over 
a background of local superfluid order does appear to lead 
to a satisfactory description of the data^ . These exper- 
iments suggest a description of the underdoped state in 
terms of a fluctuating superconducting order parameter, 
^'sc- Numerous such theories^ have been considered in 
the literature, under guises such as "preformed pairs" and 
"phase fluctuations" . Intriguing new evidence in support 
of such a "bosonic Cooper pair" approach has been pre- 
sented in recent work by Kapitulnik and collaborators^. 

On the other hand, a somewhat different view of the 
underdoped state has emerged from a second class of 
experiments. Scanning tunnelling microscopy (STM) 
studies^"^ of Bi2Sr2CaCu208+5 and Ca2_a;Naa;Cu02Cl2 
show periodic modulations in the electronic density of 
states, suggesting the appearance of a state with den- 
sity wave order. Such density wave order can gap out 
portions of the Fermi surface, and so could be responsi- 
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blc for onset of a spin-gap observed in the underdoped 
regime. However, the strength of the density wave mod- 
ulations is quite weak, and so it is implausible that they 
produce the large spin gap. The precise microscopic na- 
ture of the density wave modulations also remains un- 
clear, and they could represent spatial variations in the 
local charge density, spin exchange energy, and pairing 
or hopping amplitudes. Using symmetry considerations 
alone, all of these modulations arc equivalent, because 
they are associated with observables which are invari- 
ant under time-reversal and spin rotations. It is useful, 
therefore, to discuss a generic density wave order, and to 
represent it by density wave order parameters /?q which 
determine the modulations in the 'density' by 

Spir)=J2PQ^"^"' (1-1) 
Q 

where Q 7^ extends over some set of wavevectors. 

We will argue here that the apparent conflict between 
these two classes of experiments is neatly resolved by 
our approach. The strong vortex fluctuations suggest we 
examine a theory of superfluidity in the vicinity of a Mott 
transition. We will show that in a quantum theory of such 
a transition, which carefully accounts for the complex 
Berry phase terms, fluctuations of density wave order 
emerge naturally, and can become visible after pinning 
by impurities. 

Within the conventional Landau-Ginzburg- Wilson 
(LGW) theory'''^", with the order parameters and 
pQ at hand, a natural next step is to couple them to 
each other. Here, one begins by writing down the most 
general effective free energy consistent with the under- 
lying symmetries. A key point is that \E'sc and pq have 
non-trivial transformations under entirely distinct sym- 
metries. The density wave order transforms under space 
group operations e.g. T^, translation by the lattice vector 
a 

Ta : PQ ^ pqe*'^ " ; ^ *sc, (1.2) 

while \l/sc transforms under the electromagnetic gauge 
transformation, G, 

G-.pQ^PQ ; ^sc^^sce'^. (1.3) 

We could combine the two competing orders in a common 
'superspin' order parameter (as in the 'S0(5)' theory^^), 
but neither of the symmetry operations above rotate be- 
tween the two orders, and so there is no a priori moti- 
vation to use such a language. The transformations in 
Eqs. (1.2) and (1.3) place important constraints on the 
LGW functional of ^'sc and pQ, which can be written as 

T = Ti + [pq] + [^sc, Pq] ■ (1.4) 

Here is an arbitrary functional of \E'sc containing only 
terms which are invariant under Eq. (1.3), while ^2 is a 

function of the pq invariant under space group operations 
like Eq. (1.2). These order parameters are coupled by J^3, 



and the symmetry considerations above imply that 
only contains terms which are products of terms already 
contained in Ti and ^2- In other words, only the en- 
ergy densities of the two order parameters couple to each 
other, and neither order has an information on the local 
'orientation' of the other order. The crucial coupling be- 
tween the orders is actually in an oscillatory Berry phase 
term, but this leaves no residual contribution in the con- 
tinuum limit implicit in the LGW analysis. To leading 
order in the order parameters, the terms in Eq. (1.4) are 

•^2[pq] = ^r\Q\\pci\^ + . . . 
Q 

^3[*«c,Pq] = 5I^'IQ|I*-I>qI' + ---- (1-5) 
Q 

The phase diagrams implied by the LGW functionals 
Eqs. (1.4) and (1.5) were determined by Liu and Fisher^, 
and are summarized here in Fig 1. The transition be- 
tween a state with only superconducting order {{^ sc) 7^ 
0' (Pq) = 0) ^"^^ ^ state with only density wave or- 
der ((^sc) = 0, (pq) ^ 0) can follow one of 3 possible 
routes: (a) via a direct first order transition; (b) via an 
intermediate supersolid phase with both order parame- 
ters non-zero, which can be bounded by two second or- 
der transitions; and (c) via an intermediate "disordered" 
phase with both order parameters zero, which can also be 
bounded by two second order transitions. An important 
prediction of LGW theory is that there is generically no 
direct second order transition between the flanking states 
of Fig 1 : such a situation requires fine-tuning of at least 
one additional parameter. 

While the predictions of LGW theory in Fig 1 are ap- 
propriate for classical phase transitions at nonzero tem- 
peratures, they develop a crucial shortcoming upon a 
naive extension to quantum phase transitions at zero 
temperature. Importantly, the intermediate "disordered" 
phase in Fig 1 has no clear physical interpretation. At 
the level of symmetry, LGW theory informs us that this 
intermediate phase fully preserves electromagnetic gauge 
invariance and the space group symmetry of the lattice. 
Thinking quantum mechanically, this is not sufficient to 
specify the wavefunction of such a state. We could guess 
that this "disordered" state is a metallic Fermi liquid, 
which does preserve the needed symmetries. However, 
a Fermi liquid has a large density of states of low en- 
ergy excitations associated with the Fermi surface, and 
these surely have to be accounted for near the quantum 
phase transitions to the other phases; clearly, such ex- 
citations are not included in the LGW framework. At 
a more sophisticated level, a variety of insulating states 
which preserve all symmetries of the Hamiltonian have 
been proposed in recent years"'^^"-'^^, but these have a sub- 
tle 'topological' order which is not contained in the LGW 
framework. Such topologically ordered states have global 
(and possibly local) low energy excitations and these 
have to be properly accounted for near quantum critical 
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FIG. 1: LGW phase diagrams of the competition between 
the superconducting and density wave orders. The choice 
between the cases (a), (b), and (c) is made by the sign of 
v\Q\ and by the values of the higher order couplings in (1.5). 
The thick line is a first order transition, and the thin lines are 
second order transitions. The "disordered" phase of LGW 
theory preserves all symmetries but is not the 'topologically 
ordered' phase which also preserves all symmetries: LGW 
theory does not predict the degeneracies of the latter state 
on topologically non-trivial spatial manifolds. It is believed 
there is no such "disordered" quantum state in the underlying 
quantum system at T = 0. 



points^^; again, such excitations are entirely absent in 
the LGW framework. Thus the purported existence of a 
truly featureless intermediate "disordered" state in Fig Ic 
suggests serious shortcomings of the LGW description of 
multiple order parameters in quantum systems^'^^^^. 

Conceptually, it is useful to emphasize that the LGW 
approach to the conducting-insulating transition is quite 
the reverse of the Mott picture. Indeed, the above sym- 
metry considerations do not at all refer to transport, or 
to localization of the carriers. To address such prop- 
erties, a conventional application of the LGW method 
would be to introduce the pq and '^sc order parame- 
ters by a mean-field decoupling of a quantum model of 
bosons^^ or fermions'^®. These order parameters, once 
established, appear then as perturbations to the "quasi- 
particles" in the decoupled model. If pq becomes strong 
enough, the coincident potential barriers can localize the 



carriers. Thus in the LGW approach the Mott transition 
eventually arises out of the competing order, not vice- 
versa. 

In this paper, we will present an alternative approach 
to studying the interplay between the superfluid and den- 
sity wave ordering in two spatial dimensions. Our analy- 
sis is based upon the recent work of Senthil et al}^ (and 
its precedents^'^'^'''^^), although we will present a differ- 
ent perspective. Also, the previous work^^ was focused 
on systems at half-filling, while our presentation is appli- 
cable to lattice systems at arbitrary densities. One of the 
primary benefits of our approach is that the formalism 
naturally excludes the featureless 'disordered' states that 
invariably appear in the more familiar LGW framework. 
Our theory also makes it evident that the experimental 
observations of 'phase' fluctuations and pinned density 
wave order noted earlier are closely tied to each other, 
and, in principle, allows a compact, unified description 
of these experiments and also of the pinned density wave 
order around vortices in the superconductor^^. 

We will present our results here in the context of boson 
systems: specifically, in a general model of bosons hop- 
ping on a square lattice with short-range interactions. 
However, all the results of this paper apply essentially 
unchanged to electronic systems with short-range pair- 
ing. This is clearly the case for a toy model in which 
electrons experience a strong attractive interaction in the 
s-wave (or d-wave, though this is less obvious) channel, 
causing them to form tightly bound "molecular" Cooper 
pairs. Below the pair binding energy - i.e. the gap 
to unbound electron excitations - a description of the 
physics in terms of bosonic Cooper pairs clearly applies 
within such a model. Since we discuss universal physics 
near a bosonic superfiuid to charge ordered transition, 
the theory here clearly describes such a transition from 
a short-range paired (also called the "strong pairing" 
regime/phase) superconductor to a density- wave ordered 
phase of any model, provided the electronic gap is main- 
tained throughout the critical region. As we discuss be- 
low, a key parameter determining the character of the 
theories presented here is the number /, which for the 
lattice boson models is the average density of bosons per 
site of the square lattice in the Mott insulating state. 
Clearly, for electrons on the lattice with filling 1 — S, 
i.e. hole density S relative to half-filling, the density of 
bosonic Cooper pairs is 
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(1-6) 



We note in passing that the density of carriers in the 
superfiuid proximate to the Mott insulator is allowed to 
be different from that of the Mott insulator, as elaborated 
in Section II F. 

These features will be explicitly demonstrated in a 
companion paper^", hereafter referred to as II. We argue 
in II that a convenient phenomenological model for cap- 
turing the spin S = sector of electrons near a superfluid- 
insulator transition is the doped quantum dimer model 



4 



of the cuprate superconductors proposed by Fradkin and 
Kivelson^^. At zero hole concentration, the quantum 
dimer model^^ generically has a Mott insulating ground 
state with valence bond solid (VBS) order^^'^'^ (this or- 
der constitutes a particular realization of the general- 
ized density wave order discussed above). Above some 
finite hole density, the dimer model has a superconduct- 
ing ground state in which neutral observables preserve 
all lattice symmetries. (As we will discuss in II, the elec- 
tronic pairing in this superconductor could have a d- or s- 
(or other) wave character, depending upon microscopic 
details not contained in the dimer model). The evolu- 
tion of the phase diagram between the insulating VBS 
state and the superconductor can occur via many possi- 
ble sequences of transitions, all of which will be shown 
to be equivalent to those in the simpler boson models 
discussed in the present paper. A duality analysis of the 
doped dimer model presented in II directly gives the ex- 
pected relation, Eq. (1.6). The paper II will also extend 
the dimer model to include fcrmionic S = 1/2 excita- 
tions: this theory will be shown to have a close connection 
to other^^ U(l) or SU(2) gauge-theoretic "slave-particle" 
approaches to cuprate physics. Alternative formulations 
appropriate to long-range paired superconducting states 
(i.e. those with gapless nodal quasiparticle excitations) 
will also be discussed in II. 

Let us turn, then, to a model of bosons hopping on 
a square lattice, with short-range interactions,^^ and at 
a density / per site of the square lattice. Our primary 
technical tool will be a dual description of this boson 
system using vortex degrees of freedom^'''"^''. Briefly, the 
world lines of the bosons in three-dimensional spacetime 
are reinterpreted as the trajectories of vortices in a dual, 
classical, three-dimensional 'superconductor' with a dual 
'magnetic' field oriented along the 'time' direction. The 
dual superconducting order parameter, tp, is the creation 
and annihilation operator for vortices in the original bo- 
son variables. The density of the vortices in ip should 
equal the density of the bosons, and hence the dual 'mag- 
netic' field acting on the tjj has a strength of / flux quanta 
per unit cell of the dual square lattice. 

It is valuable to characterize the phases in Fig 1 in 
terms of the dual field ip- It is easy to determine the 
presence/absence of superfluid order (which is associated 
with the presence/absence of a condensate in '^sc)- This 
has a dual relationship to the dual 'superconducting' or- 
der and is therefore linked to the absence/presence of a 
condensate in ^: 

Superfluid : (ijj) = 
Insulator : (V-) 7^ 0. (1.7) 

A condensate in tjj implies a proliferation of vortices in 
^sc, and hence the loss of superfluid order. Strictly 

speaking, the loss of superfluidity does not necessarily 
require the appearance of a condensate of elementary 
vortices. It is sufficient that at least a composite of n 
vortices condense, with (V'") 7^ 0, where n is a positive 
integer. The insulating states with n > 1 and {tp) = 



have topological order . For simplicity, we will not con- 
sider such insulating states here, although it is not diffi- 
cult to extend our formalism to include such states and 
the associated multi- vortex condensates. 

The characterization of the density wave order in terms 
of the vortex field, ip, is more subtle, and a key ingredient 
in our analysis. We have noted that the field ip experi- 
ences a dual magnetic field of strength / flux quanta per 
unit cell. Consequently, just as is familiar from the Hofs- 
tadter problem of electrons moving in a crystal lattice in 
the presence of a magnetic field^*'^^, the ip Hamiltonian 
transforms under a projective representation of the space 
group''^'''^ or a PSG. A central defining property of a 
projective representation is that the operations T^, Ty of 
translation by one lattice site along the x, y directions do 
not commute (as they do in any faithful representation 
of the space group), but instead obey 

TxTy = uTyTx, (1-8) 

where 

oj = e^"*^ . (1.9) 

The PSG also contains elements corresponding to all 
other members of the square lattice space group. Among 
these are R^^f (rotation by angle tt/2 about a site of the 
dual lattice i.e. the sites of the lattice upon which the 
vortex field tjj resides) and Z^"*^', /^"^' (reflections about 
the X, y axes of the dual lattice); these obey^^ group mul- 
tiplication laws identical to those in the ordinary space 
group, e.g. 

rp pdual pdualrri— 1 

All the elements of the PSG are obtained by taking ar- 
bitrary products of the elements already specified. 

A useful and explicit description of this projective rep- 
resentation is obtained by focusing on the low energy and 
long wavelength fluctuations of ip at boson density 

/ = ^ (1-11) 

where p and q are relatively prime integers. We assume 
that the average density in the Mott insulator has the 
precise commensurate value in Eq. (1.11). In most of 
the paper we will also assume that the superfluid has the 
density value in Eq. (1.11), but it is not difficult to ex- 
tend our analysis to allow the density of the superfluid 
(or supersolid) to stray from commensurate values to ar- 
bitrary incommensurate values: this will be discussed in 
Section II F. As we will review in Section II, the spectrum 
of Ip fluctuations in the phase with {ip) = has q distinct 
minima in the magnetic Brillouin zone. We can use these 
minima to define q complex fields, i? = 0, 1, . . . , q — 1 
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which control the low energy physics. These q vortex 
fields play a central role in all our analyses, as they allow 
us to efBciently characterize the presence or absence of 
superfluid and/or density wave order. In a convenient 
Landau gauge, the q vortex fields can be chosen to trans- 
form under and Ty as 



(1.12) 



Here, and henceforth, the arithmetic of all indices of the 
ipi fields is carried out modulo e.g. ipq = <po- The 
action of all the PSG elements on the ipe will be specified 
in Section II B; here we also note the important transfor- 
mation under i?*^*)?' 

7r / Z 
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(1.13) 



which is a Fourier transform in the space of the q fields. 
It is instructive to verify that Eqs. (1.12) and (1.13) obey 
Eqs. (1.8) and (1.10). 

The existence of q degenerate species of vortices in the 
superfluid is crucial to all the considerations of this pa- 
per. At first sight, one might imagine that quantum tun- 
nelling between the different vortex species should lift the 
degeneracy, and we should only focus on the lowest en- 
ergy linear combination of vortices. However, while such 
quantum tunnelling may lead to a change in the preferred 
basis, it does not lift the (/-fold degeneracy. The degener- 
acy is imposed by the structure of the PSG, in particular 
the relation Eq. (1.8), and the fact that the minimum 
dimension of a PSG representation is q. 

Using the transformations in Eqs. (1.12) and (1.13), 
we can now easily construct the required density wave 
order parameters pq as the most general bilinear, gauge- 
invariant combinations of the ip^ with the appropriate 
transformation properties under the square lattice space 
group. The density wave order parameters appear only 
at the wavevectors 



Qmn = 27r/ (m,n), 
where m, n are integers, and take the form 



(1.14) 
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p„„ = pQ_ = 5(|Q„„|)u;™"/2^^|^^+„u;^'". (1.15) 

Here S{Q) is a general 'form-factor' which cannot be 
determined from symmetry considerations, and has a 
smooth Q dependence determined by microscopic details 
and the precise definition of the density operator. It is 
easy to verify that — p-m,-n, and from Eqs. (1-12) 
and (1.13) that the space group operations act on 
just as expected for a density wave order parameter 
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(1.16) 



So by studying the bilinear combinations in Eq. (1.15), 
we can easily characterize the density wave order in terms 
of the Iff. It is worth reiterating here that the pmn or- 
der parameters in Eq. (1.15) are generic density wave 
operators, and could, in principle, be adapted to obtain 
the actual density of either the bosons or the vortices. 
Indeed, as we noted above Eq. (1.1), we can obtain in- 
formation on arbitrary observables invariant under spin 
rotations and time reversal, including the local density of 
states. The different specific observables will differ only 
in their values of the form factor S{Q), which we do not 
specify in the present paper. Note also that because the 
^Tr/f rotation is about a dual lattice site, the last relation 
in Eq. (1.16) implies that the Fourier components pq are 
defined by taking the origin of the spatial co-ordinates on 
a dual lattice site. Upon performing the inverse Fourier 
transform from the pmn to real space (as in Eq. (2.26) 
below) the resulting density gives a mcasTire (up to a fac- 
tor of the unknown S{Q)) of the vortex density on the 
dual lattice sites, and of the boson density on the direct 
lattice sites. 

For completeness, we note that the characterizations 
in Eq. (1.7) of the superfluid order in terms of the vortex 
fields has an obvious expression in terms of the (ff. 

Superfluid : (ipe) = for every £ 
Insulator : (tp^) ^ for at least one £. (1.17) 

The key relations in Eqs. (1.15) and (1.17) allow us 
to obtain an essentially complete characterization of all 
the phases of the boson model of interest here: this is a 
major reason for expressing the theory in terms of the 
vortex fields ip^. Furthermore, the symmetry relations 
in Eqs. (1.12) and (1.13) impose strong constraints on 
the effective theory for the (fi fields whose consequences 
we will explore in Section II C. Indeed, even before em- 
barking upon a study of this effective theory, we can al- 
ready see that our present theory overcomes one of the 
key shortcomings of the LGW approach. We note from 
Eq. (1.17) that the non-superfluid phase has a condensate 
of one or more of the (pi. Inserting these condensates into 
Eq. (1.15) we observe that (pmn) 7^ for at least one non- 
zero Qm„ i-P'- the non-superfluid insulator has density 
wave order. Consequently the intermediate "disordered" 
phase in Fig Ic has been automatically excluded from the 
present theory. Instead, if the (p£ condensate appears in 
a second-order transition (something allowed by the gen- 
eral structure of the theory), we have the possibility of 
a generic second-order transition from a superfluid phase 
to an insulator with density wave order, a situation which 
was forbidden by the LGW theory.'''' We summarize the 
possible phase diagrams of the (p^ theory in Fig 2. 

In addition to characterizing modulations in the "den- 
sity" in terms of the p„m , the (p£ flclds can also determine 
modulations in the local vorticity. This is a measure of 
circulating boson currents around plaquettes of the di- 
rect lattice^"', and has zero average in all phases because 
time-reversal invariance is preserved. However, in the 
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FIG. 2: Possible phase diagrams in the vortex theory of the 
interplay of superfluid and density wave order. Conventions 
are as in Fig 1. Notice that the objectionable "disordered" 
phase of the LGW theory has disappeared, and is replaced 
by a direct second order transition. Also, the first order tran- 
sition in (a) is "fluctuation induced", i.e. unlike the LGW 
theory of Fig 1, the mean field vortex theory only predicts a 
second order transition, but fiuctuations could induce a weak 
first order transition. 



presence of an applied magnetic field, this constraint no 
longer applies. The magnetic field, of course, induces a 
net circulation of boson current, and the total number 
of vortices (or anti-vortices) is non-zero. However, using 
a reasoning similar to that applied above to the density, 
the PSG of the vortices implies there are also modula- 
tions in the vorticity at the Qmn wavevectors. We can 
define a corresponding set of Vmn which are Fourier com- 
ponents of the vorticity at these wavevectors. Details of 
this analysis, and explicit expressions for the Vmn appear 
in Appendix A. 

An explicit derivation of the above vortex theory of 
the Mott transition of the superfluid appear in Section II 
where we consider a model of bosons on the square lat- 
tice. Here we derive an effective action for the q vortex 
species in Eqs. (2.19) and (2.21) which is invariant un- 
der the PSG, and which controls the phases and phase 
diagrams outlined above. 

A complementary perspective is presented in Sec- 
tion III where we formulate the physics using the direct 



boson representation. In particular, following Senthil et 
al.^^, we explore the possibility that the LGW forbidden 
superfluid-insulator transitions discussed above are asso- 
ciated with the boson fractionalization. The structure of 
the theory naturally suggests a fractionalization of each 
boson into q components, each with boson number l/q. 
However, the PSG places strong restrictions on the frac- 
tionalized boson theory. On the rectangular lattice, such 
a theory is consistent for all q but only for certain insulat- 
ing phases obtained for restricted parameter values. On 
the square lattice, with the additional W^^^^ element of 
the PSG, a fractionalized boson theory is permitted only 
for a limited set of q values. Specifically, it is required 
that the ipe vortex fields transform under a permutative 
representation of the PSG. We define a permutative rep- 
resentation as one in which all group elements can be 
written as AP, where A is a unitary diagonal matrix 
(i.e. a matrix whose only non-zero elements are complex 
numbers of unit magnitude on the diagonal), and P is 
a permutation matrix. For the representation defined in 
Eqs. (1.12) and (1.13), and Ty are of this form, while 
R'^^jf is not. When a permutative representation exists, 
it is possible to globally unitarily transform the ife fields 
to a new basis of Q fields, which realize the permutative 
representation. For some range of parameters - values of 
the quantum tunnelling terms between the vortex flavors 
- one can see that the PSG is expected to be elevated at 
the critical point to include an emergent C/(l)' (modulo 
the global U{1) gauge invariance) symmetry of indepen- 
dent rotations of the Q fields. This emergent symmetry 
is the hallmark of fractionalization, and corresponds to 
emergence of conserved U(l) gauge fluxes in the direct 
representation. Such an emergent symmetry is clearly 
not generic, and this turns out to place quite a restrictive 
condition on the existence of a fractionalized l/q boson 
representation, as we discuss in detail in Section HI, and 
in a number of appendices. 

Section IV applies our theory to the STM observations 
of local density of states modulations in zero and non- 
zero magnetic field. A virtue of our approach is that it 
naturally connects density wave order with vorticity and 
allows a unified description of the experiments in zero 
and non-zero magnetic field. As we will discuss, there 
are several intriguing consequences of our approach, and 
we will analyze prospects for observing these in future 
STM experiments in Sections IV and V. 

We briefly mention other recent works which have ad- 
dressed related issues. Zaanen et al.^'^ have examined the 
superfluid to insulator transition from a complementary 
perspective: they focus on the dislocation defects of a 
particular insulating solid in the continuum, in contrast 
to our focus here on the vortex defects of the superfluid. 
We do consider the 'melting' of defects in the solid in Sec- 
tion HI C, but the underlying lattice plays a crucial role in 
our considerations. In a work which appeared while our 
analysis was substantially complete, Tesanovic"^^ has ap- 
plied the boson-vortex duality to Cooper pairs and con- 
sidered the properties of vortices in a dual 'magnetic' 
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field. 



II. 



DUAL VORTEX THEORY OF BOSONS ON 
THE SQUARE LATTICE 



We consider an ordinary single-species boson model on 
the square lattice. The bosons are represented by rotor 
operators (pi and conjugate number operators hi where 
i runs over the sites of the direct square lattice. These 
operators obey the commutation relation 



(2.1) 



A simple boson Hamiltonian in the class of interest has 

the structure 

n = -t^ cos (^Aa^i - 2irgia^ + ^ V(ni) 



^^Aijfiihj + . . 



(2.2) 



where the interaction V{'h) has the on-site terms V{h) — 
—fih + Uh{h — l)/2. The Aij are repulsive off-site in- 
teractions, and can also include the long-range Coulomb 
interaction. We will focus on the case of short-range Aij , 
and note the minor modifications necessary for the long- 
range case. The general structure of our theory also per- 
mits a variety of exchange and ring-exchange terms, such 
as those in the studies of Sandvik et al.^^. These off-site 
or ring exchange couplings are essential for our analysis, 
as they are needed to stabilize insulating phases of the 
bosons away from integer filling. Nevertheless, many as- 
pects of our results are independent upon the particular 
form of these couplings; their specific form will only infiu- 
ence the numerical values of the non-linear couplings that 
appear in our phenomenological ac;tions (such as those in 
Eq. (2.21)). The index a extends over the spatial direc- 
tions X, y, while we will use indices /x, v, A to extend 
over all three spacctimc directions .t, y, r. The symbol 
Aa is a discrete lattice derivative along the a direction: 
Aa(f>i = (f)i+a — (t> (and similarly for A^). We have also in- 
cluded a static external magnetic field represented by the 
vector potential gia for convenience. This is a uniform 
field which obeys 



ivxA^gix = hSfj 



(2.3) 



where h is the strength of the physical magnetic field 
(which should be distinguished from the dual "magnetic" 
flux / discussed in Section I). 



A. Dual lattice representation 

We proceed with a standard duality mapping, follow- 
ing the methods of Refs. 25-27 and the notational con- 
ventions of Ref. 37. We represent the partition function 



as Fcynman integral over states at a large number of in- 
termediate time slices, separated by the interval At. The 
intermediate states use a basis of hi and 4>i at alternate 
times. The hopping term in Ti. acts between (j)i eigen- 
states, and we evaluate its matrix elements by using the 
Villain representation 



exp ^fArcos ^Aq^j 



E 



exp 



2tAT 



+ iJiaAa4>i - 2mJiagia ] ■ (2-4) 



We have dropped an unimportant overall normalization 
constant, and will do so below without comment. The 
Jia are integer variables residing on the links of the direct 
lattice, representing the current of the bosons. 

After integrating over the (f)i on all sites and at all 
intermediate times, the partition function becomes 




2e2 E ^"^^t" 



Ar^AijJirJjT - ^wigi^Jifj, \ 



(2.5) 



where Jj^ = (n^, Jix, Jiy) is the integer- valued boson cur- 
rent in spacetime, i now extends over the sites of the cu- 
bic lattice, and we have chosen At so that = tAr = 
1/UAt, and H = fi/U + 1/2. We now solve the con- 
straint in Eq. (2.5) by writing 



A„A 



aX 



(2.6) 



where a labels sites on the dual lattice, and Aax is an 
integer- valued gauge field on the links of the dual lattice. 
We promote Aa/^ from an integer-valued field to a real 
field by the Poisson summation method, while "soften- 
ing" the integer constraint with a fugacity yy. It is con- 
venient to make the gauge invariance of the dual theory 
explicit by introducing an angular field da on the sites of 
the dual lattice, and mapping 2TTAa^ — > ^irAa^i — A^i?a- 
The operator e*''" is then the creation operator for a vor- 
tex in the boson phase variable (jji. These transformations 
yield the dual partition function 



JJ j dAau, j dt?aexp|^ 

-ihY,{^r^a-2T:Aar)\. (2.7) 
a / 
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We have not explicitly displayed the Ay term above, and 
assumed it has been absorbed into a renormalized value 
of e^.^^ This is the theory of a dual vortex boson repre- 
sented by the angular rotor variable i?a) coupled to a dual 
gauge field Aa^. The Berry phase term proportional to 
h is precisely the constraint that the rotor number vari- 
able (canonically conjugate to i?o) takes values which are 
integers plus h (see Ref. 27): so there is a background 
density of vortices of h per site which is induced by the 
magnetic field acting on the direct bosons. 

The remainder of Section II will consider only the case 
h = 0; the analog of the h case will appear later in 
the dimer model analyses of II. 

We obtain the dual theory in its final form on the cubic 
lattice by replacing the field e*''° by a "soft-spin" dual 
vortex field tpa, which yields 

a \ 



afj. 



(2.8) 



The last line in Eq. (2.8) is the effective potential for 
the complex vortex field ipa- Increasing the parameter s 
scans the system from the insulating solid at s <C to 
the superfiuid at s » i.e. the system moves from right 
to left in the phase diagrams of Fig 2. 



B. Symmetries 

We will now present a careful analysis of the symme- 
tries of Eq. (2.8), with the aim of deducing general con- 
straints that must be obeyed by the low energy theory 
near the superfluid-to-insulator transition. 

We begin in the superfiuid regime with ,s large, so that 
(V'a) = 0. The direct boson density is the r component of 
the dual 'magnetic' fiux e^v\A^A\, and for s large in the 
action in Eq. (2.8) it is clear that the saddle point of the 
A^ fluctuations occurs at Aa^ 



with eui,\Ai,Aax 



HSfiT- We want this boson density to be the value / in 
Eq. (1.11), and so we should choose H = f. 

We now wish to examine the structure of ipa fluctua- 
tions about this saddle point. Let us choose the Landau 
gauge with Aar = Aax = and 



(2.9) 



here is the x co-ordinate of the dual lattice a. We now 
need to determine the ipa spectrum in a background Aa^ 
held. 



dual 



The basic symmetry operations arc T^, Ty, and ' 
introduced in Section I. The action of these operators on 



tpa = tp{ax,ay) required to keep the Hamiltonian invari- 
ant is 



Ty 

Tx 
odual 

-^77/2 



il){ax,ay) ip{ax,ay - 1) 

ipi^ax, ay) %l>{ax - 1, ay)ijj°'y 

^p{ax, tty) ^p{ay, -a,;)u;"»=°'' (2.10) 



Notice that Eqs. (1.8) and (1.10) are obeyed by the above. 

It is also useful to collect the representation of 
Eq. (2.10) in momentum space. Implying that all mo- 
menta are reduced back to the extended Brillouin zone 
with momenta —7r<kx,ky<7rwe flnd 



R. 



.dual 
•7r/2 



ip{kx,ky) ip{kx,ky)e~'''y 
ij:{kx, ky) V(fcx, ky - 27r/)e-''=» 
i^{kx,ky)^ (2.11) 

1 

- i;{ky + 2TTnf,-kx-2TTmf)oj-"''' 

^ m,n=0 



The most important tpa fluctuations will be at mo- 
menta at which the spectrum has minima. It is not diffi- 
cult to show from the above symmetry relations that any 
such minimum is at least gr-fold degenerate: Let |A) be a 
state at a minimum of the spectrum. Because, the oper- 
ator Ty commutes with the Hamiltonian, this state can 
always be chosen to be an eigenstate of Ty, with eigen- 
value e"'*^" . Now the above relations imply immediately 
that Tx\A) is also an eigenstate of the Hamiltonian with 
Ty eigenvalue e~'^'^yu)~^. Also, because the Ty eigenvalue 
is distinct from that of |A), this state is an orthogonal 
eigenstate of the Hamiltonian. By repeated application of 
this argument, we obtain q orthogonal eigenstates of the 
Hamiltonian whose Ty eigenvalues are e~^''y times integer 
powers oiu}~^. Because we are working with the gauge of 
Eq. (2.9), our Hamiltonian couples momenta {kx,ky) to 
momenta {kx ± 27r/, ky). It is therefore advantageous to 
look at the spectrum in the reduced Brillouin zone with 
—ir/q < kx < ir/q and — tt < ky < tt. For the nearest- 
neighbor model under consideration here, the minima of 
the spectrum are then at the q wavevectors {0,2n£p/q) 
with £ = 0, . . .q—1. Let us label the eigenmodes at these 
wavevectors tpg. We therefore have to write down the 
field theory in terms of these q complex fields ipe. 

It is useful, especially when analyzing the influence of 
to make the above symmetry considerations ex- 



t/2 ■ 

plicit. In the extended Brillouin zone, for the nearest- 
neighbor model under consideration here, let us label the 
q'^ states at the wavevectors {2'jTmp/q, 2imp/q) by |m, n). 
Then one of the minima of the spectrum corresponds to 
the state 



9-1 

bo) = Y '^'"I'^'O) 

TO=0 



(2.12) 



where the Cm are some complex numbers. Then, by oper- 
ation of Tx on |(po) we obtain the q degenerate eigenstates 
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as 



9-1 



(2.13) 



m=0 



Now let us consider the action of states 
in Eq. (2.13). Using Eq. (2.11), and after some simple 
changes of variables we obtain 



odual 

•^^77/2 



(2.14) 

Now, because commutes with H, the right-hand- 

side of Eq. (2.14) must be a linear combination of the \Lpi) 
states in Eq. (2.13). The matrix elements of the rota- 
tion operator are then given by {^e\R'^i2\'-Pi) = clo~^^ , 
where c = | cJ^/Cm'^'""* is independent of i and i'. 
For the nearest-neighbor model under consideration here, 
it can be easily checked that the c^'s are invariant under 
a Fourier transform such that c = l/y^- Hence we find 



1 



V ^ i'=0 



(2.15) 



Our discussion above has now established that the low 
energy vortex fields must have an action invariant under 
the transformations in Eqs. (1.12) and (1.13). In a sim- 
ilar manner we can also determine the transformations 
associated with the remaining elements of the square lat- 
tice space group. These involve the operations I^"**' and 
7^"*' which are reflections about the x and y axes of the 
dual lattice. Under these operations we find 



rdual 



7-dual 



* 



(2.16) 



Finally it is interesting to consider the point inversion 
operator = (i?^^^')^, with 

rdual , ,^ 

As in the ordinary space group we have 



rdual 



dual rdual 

y 



rdual rdual 
y X 



(2.17) 



(2.18) 



will be represented by the continuum non-compact U(l) 
gauge field a^p/(27r), where a is the lattice spacing. 

First, we consider quadratic order terms about the sad- 
dle point of Eq. (2.8). The most general action has the 
familiar terms of scalar electrodynamics 

So = f cfrdT (j2 [K^/^ - i^M'' + 



(2.19) 



We have rescaled the coupling e here by a factor of 27r 
from Eq. (2.8). 

Next, we consider terms which are quartic in the ipi, 
but which contain no spatial or temporal derivatives. 
These will be contained in the action Si. We discuss 
two approaches to obtaining the most general quartic in- 
variants. The first is the most physically transparent, 
but turns out to be eventually inconvenient for explicit 
computations. In this approach we use density operators 
defined in Eq. (1.15), and their simple transformation 
properties in Eq. (1.16), to build up quartic invariants. 
In particular, the quartic invariants are only quadratic 
in the Pmm and we need only the most general quadratic 
term invariant under Eq. (1.16). This has the form 



„ /q/2 n 

Si = d^^dr EE^^^O 

\n=Om=0 



1 2 I 1 2 

Pnm\ ~l~ ml 



(2.20) 



However, not all the invariants above are independent, 
and there are often linear relations between them - this 
reduces the number of independent coupling constants 
Xnm- Determining the linear relations between the cou- 
plings turns out to be inconvenient, and we found it easier 
to proceed by the second method described below. 

In the second approach, we first impose only the con- 
straints imposed by the translation operations Tx, Ty. By 
inspection, it is easy to see that the most general quartic 
term invariant under these operations has the structure 



Si 



(2.21) 



C. Continuum field theories 

We have established that fluctuations of the vortex 
fields about the saddle point of Eq. (2.8) in Eq. (2.9) 
transform under a projective representation of the square 
lattice space group which is defined by Eqs. (1.12), (1.13), 
and (2.16). In this section we will write down the most 
general continuum theory of the fields which is in- 
variant under these projective transformations. The ac- 
tion should include fluctuations in Aa^ about A^n - this 



Here the integers i,jn,n,... range implicitly from to 
q—1 and all additions over these integers are taken mod- 
ulo q. Imposing in addition the reflection operations in 
Eq. (2.16), and accounting for the internal symmetries in 
Eq. (2.21), it is easy to show that the couplings 7„,„ can 
always be taken to be real and to obey the relations 



^mn ~ Tm— 2n, — n* 



(2.22) 
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We have so far not yet imposed the constraints impHed 
by the lattice rotation R'^^f- Consequently, Eqs. (2.21) 

and (2.22) define the most general quartic terms for a 

system with rectangular symmetry, which may be appro- 
priate in some physical situations. By explicit solution of 
the constraints defined by Eq. (2.22), we found that the 
quartic terms are dc^tc^mined by Afreet independent real 
coupling constants, where 

AT^ect = ^^+^)^(^ + ^) for q = 2n, 2n + 1, (2.23) 

with n a positive integer. Note that the number of cou- 
plings grows quite rapidly with increasing q: Afreet ~ 

Finally, let us consider the consequences of the 
symmetry. By applying Eq. (1.13) to Eq. (2.21) we find 
that <Si remains invariant provided the ^mn obey the re- 
lations 



q ^ 

mn 



^ — [n{m—n)-\-n{m—n)] 



(2.24) 



Naively, there arc relations implied by Eq. (2.24), but 
they are not all independent of each other. By explic- 
itly solving the relations in Eq. (2.24) for a range of q 
values we found that the number of independent quartic 
coupling constants for a system with full square lattice 
symmetry is 



_( (n-l-1)^ for = 4n, 4n -I- 1 
~ \{n+ l)(n + 2) for q = in + 2, An + 3 



(2.25) 

with n a positive integer. The additional restrictions of 

square relative to rectangular symmetry reduce the num- 
ber of independent coupling constants by roughly half at 
large q: iVgquare ^ 9^/16. 



D. Mean field theory 

This section will examine the mean-field phase dia- 
grams of the general theory <So -|- Si proposed in Sec- 
tion lie. From the discussion in Section I, it is clear 
that such a procedure can yield a direct second order 
transition from a superfluid state (with {(pi) — 0) to an 
insulating state with density wave order (with (tp^) ^ 0). 
We are interested in determining the possible configura- 
tions of values of the (pe and the associated patterns of 
density wave order for a range of q values (in this sub- 
section, we will simply write {(pe) as (pe because there is 
no distinction between the two quantities in mean field 
theory). After determining the ipe, we determined the 
Pmn by Eq. (1.15) (using a Lorentzian for the form fac- 
tor S{Q) = 1/(1-1- Q^)), and then computed the density 
wave order using 



g-l 



(2.26) 



We evaluated Eq. (2.26) at r values corresponding to 
the sites, bonds, and plaquettcs of the direct lattice, 
and plotted the results in the square lattice figures that 
appear below. In the formalism of Section II, r values 
with integer co-ordinates correspond to the sites of the 
dual lattice, and hence to plaqucttes of the direct lattice. 
The value of Sp{r) on such plaquette co-ordinates can be 
considered a measure of the ring-exchange amplitude of 
bosons around the plaquette. By a similar reasoning, r 
values with half-odd-integer co-ordinates represent sites 
of the direct lattice, and the values of 6p{r) on such sites 
measure the boson density on these sites. Finally, r val- 
ues with r^ integer and ry half-odd-intcger correspond to 
horizontal links of the square lattice (and vice versa for 
vertical links), and the values of 6p{r) on the links is a 
measure of the mean boson kinetic energy; if the bosons 
represent a spin system, this is a measure of the spin 
exchange energy. 

Our results appear for a range of q values in the fol- 
lowing subsections. 



1. q = 2 

After imposing the symmetry relations in Eqs. (2.22) 
and (2.24), the quartic potential in Eq. (2.21) (defined 
by iSi = / (frdTCi) is 



700 



A = ^ (|<PoP + Wi?Y + ^ (Vo^t - vl'Pif (2.27) 
We now make the change of variables 

Co + Ci 



Vo 



■ 



V2 

■ Co-Ci 
— 
^2 



(2.28) 



As we will discuss in detail in Section III A, the new vari- 
ables have the advantage of realizing a permutative rep- 
resentation of the PSG, i.e. all representation matrices 
of the PSG have only one non-zero unimodular element 
in every row and column. From Eqs. (1.12), (1.13), and 
(2.28) it is easy to deduce that in the Q variables. 



-i 

1 



odual 
^7r/2 



1 

1 

n p-*7r/4 

W4 ) • (2-29) 



The action in Eq. (2.27) reduces to 



A = ^ (ICoP + ICiH' - ^ (ICol^ - ICln^ (2.30) 

The result in Eq. (2.30) is identical to that found in ear- 
lier studies^^'^* of the q = 2 case. 

Minimizing the action implied by Eq. (2.30), it is evi- 
dent that for 701 < there is a one parameter family of 




FIG. 3: Mean field phase diagram of So + Si defined in 
Eqs. (2.19) and (2.21). The couplings obey Eqs. (2.22) and 
(2.24). We assumed s = —1 to obtain an insulating state, 
and 700 = 1. Plotted are the values of Sp{r), defined as in 
Eq. (2.26), on the sites, links and plaquettes of the direct 
lattice: the lines represent links of the direct lattice. As dis- 
cussed in the text, these values represent the boson density, ki- 
netic energy, and ring-exchange amplitudes respectively. The 
saddle point values of the fields in the states are shown in 
Eq. (2.31). The choice between the states B and C is made 
by a eighth order term in the action. 



gauge- invariant solutions in which the relative phase of Co 
and Ci remains undetermined. Computing the values of 
Sp{r) establishes that this phase is physically significant, 
because it leads to distinguishable density wave modula- 
tions. This indicates that the phase will be determined by 
higher order terms in the effective potential. As shown in 
earlier work^^, one needs to go to a term which is eighth 
order in the ipe before the value of arg(Co/Ci) is pinned 
at specific values. In the present formalism, the required 
eighth order term is Pq iPi q, which is invariant under all 
the transformations in Eq. (1.16); this contains a term 

~(Cocr)'+c.c. 

Our specific results from the mean field theory for 
q — 2 are summarized in Fig 3. The state (A) has an or- 
dinary charge density wave (CDW) at wavevector (7r,7r). 
The other two states are VBS states in which all the sites 
of the direct lattice remain equivalent, and the VBS or- 
der appears in the (B) columnar dimer or (C) plaquette 
pattern. The saddle point values of the fields associated 
with these states are: 

(A) : Co 7^ , Ci = or Co = , Ci 7^ 0. 

(B) : Co = e'""/'Ci 0. 

(C) : Co = e'("+i/2)-/2Ci ^ 0, (2.31) 

where n is any integer. At quartic order, the states B 
and C are degenerate with all states with Co = e^^'Pi, 
with d arbitrary; the value of 6 is selected only by the 
eighth order term. 
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FIG. 4: As in Fig 3 but for g = 3 

2. (7 = 3 

Now there are 3 ipi fields, and the quartic potential in 
Eq. (2.27) is replaced by 

A = ^(I^0P+|^1P + |^2P)' 

I 701 /'**2i **2i **2i 

+ — (Vo'^i'/'2 + 'y='i'y52'^o + V2Vo'^i +C-C. 

- 2\^o\'\^i\' - 2|^in^2p - 2|^2n¥'oP)(2.32) 

We show in Appendix E that there is now no transfor- 
mation of the ip£ variables which realize a permutative 
representation of the PSG, and in which Eq. (2.32) may 
take a more transparent form. The results of the mini- 
mization of Sq + Si for the insulating phases are shown 
in Fig 4. The states have stripe order, one along the 
diagonals, and the other along the principle axes of the 
square lattice. Both states are 6-fold degenerate, and the 
characteristic saddle point values of the fields are 

(A) : (fio , ipi ^ (p2 = or 

(B) : ^0 - '/'I - e±2W3<^2 

and permutations. (2.33) 

3. g = 4 

For q = 4, there are 4 independent quartic invariants in 
Eq. (2.21), and they can be combined into the following 
forms. 

/1=P0,0 = (I^0P + |^1P + |^2P + |^3P)' 

h = l^ol'lV'll' + \Vl\'\^2\^ + \V2f\^3\' + ImflVof 

-fofif^fa - vlvWaVQ - flvlvo'fi - <P3'pS¥'i<P2 

I4 = 1^0^*2 {fl + 'fl) + C-C- + {^0 + Vl) + C-C. 

-4|^o|'|^2|'-4|v:>i|V3l' (2.34) 

Now as for q = 2, the invariants for g = 4 do simplify by 
redefinitions of the fields analogous to Eq. (2.28), which 
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realize a permutative representation of the PSG. We de- 
fine 



<^3 



(Co + Cl)/^^ 
(Co " Ci)/V2 

(C3 + C2)/%/2 
(C3 - C2)/V2. 



(2.35) 



The PSG operations on the Q can then easily be shown 
to take the permutative form specified in Eq. (C5) (in 
a slightly different notation). In the new variables, the 
quartic potential can be written in the following form 



£4 = 



^(ICol 
^^i(ICo|' + 

^^2(|ConC2 



icir + ic2r + ic3i 



IC2'^^ 



(ICiP- 



,21 ) W^l 
2 I lA |2|/- |2 



ic3n 



ICillCsl^ 



*2/-2 
2 



C2^C? 



+ C2'C3' + C3'C| + C3'Co' + Co*'C|) (2.36) 

where u, vi^2 and A are couphng constants. Notice that at 
A = the action is independent of the relative phases of 
the Cij and so there is a 3 parameter family of degenerate 
states. It is therefore convenient to first determine the 
minima at A = (which can be done analytically), and 
to then determine the fate of the minima so found at 
small A. The results of such a procedure are summarized 
in Fig 5. The states in Fig 5 are characterized by the 
following parameter values 



{A and B) 

C 
D 

E 



One of Co or C2 non-zero, one of Ci or Cs 

non-zero, magnitudes of non-zero Q equal 

Only one of \Q\ non-zero 

Either Co - C2 = and |Ci | = ICsl 0, 

or Ci = C3 = and |Co| = IC2I 7^ 

All \Q\ equal and non-zero (2.37) 



By inserting the saddle point values of state D in 
Eq. (2.36) we observe that this state has a continuous 
degeneracy with Ci = e'^Cs (for the first choice of the 
saddle point) and 9 arbitrary. We expect that higher or- 
der terms in the action will select the value of 9, and we 
show the higher symmetry states at special values of 9 in 
Fig 6. 

The phases in Fig 5 all survive for a finite range of 
A. However, for large enough |A| a number of additional 
phases, with distinct patterns of lattice symmetry break- 
ing, appear. We have not determined the full phase dia- 
gram in the three-dimensional vi, V2, A space, but here 
merely list the additional phases that were found in our 
numerical analysis. The new large |A| phases are shown 
in Fig 7, and the typical saddle points had the following 
characteristics: 
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FIG. 5: As in Fig 3, but for q — 4 with the quartic potential 
as in Eq. (2.36). We have set u = 1, and then determined 
the phase diagram as function of vi and V2 for infinitesimal 
values of A. The dashed lines represent the vi and V2 axes, 
while the full lines are phase boundaries. 





FIG. 6: Higher symmetry cases of the state D in Fig. 5. With 
the angle 6 defined as in the text, state Di corresponds to 
9 — n7r/2, and the state D2 to 9 = {n + l/2)7r/2 (n integer). 



F : Co = C3, Ci = C2, ICoIt^ICiI 
G : Co = *C2, |Ci|7^lCo|,arg(Co) = 

and C3 = 
H : C3 = <o, Ci = <2, IColy^lCil 
I ■■ Ci=<3, ICi|y^lC2|,arg(C2) = 

and Co = 0. 

4- General q 



ai-g(Ci 



arg(Ci) 



(2.38) 



The insulating states obtained above clearly have the 
dimensions of their unit cells constrained by the value of 
q, and we will discuss these constraints more explicitly 
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FIG. 7: Additional insulating mean-field states for q = 4 with 
the quartic potential as in Eq. (2.36). These states appear for 
larger values of |Aj and the field configurations are described 
in the text. 



here with special attention to the value q — 16. From 
Eq. (1.6) and the analysis of II, this value of q is of par- 
ticular relevance to the cuprates at the hole density of 
1/8. From Eq. (2.25), we note that the q = 16 case 
has 25 independent quartic coupling constants. The full 
parameter space of insulating states is therefore of im- 
mense complexity, and we will not attempt to map out 
any phase diagram. Instead we shall examine the spatial 
structure of a few simple configurations of the tpe- 

Let us assume that the insulating state has a unit cell 
of size a X b lattice sites (a, b integers). We now discuss 
the general constraints on the values of a and b. From 
the allowed values of the wavevectors in Eq. (1.14) we 
see that generically we will have a = q and b — q. How- 
ever, it could be that the (pe take special values so that 
some of the pmn vanish. In this case, the periods a and 
b could be any integer divisor of q. However, an inter- 
esting property of Eq. (1.15) is that it is not possible to 
simultaneously reduce the values of both a and b. The 
non-commutativity of and Ty in Eq. (1.8) effectively 
imposes an "uncertainty relation" between their values 
so that the product ab must be a multiple of q: intu- 
itively, we see this as a requirement that the unit cell 
must contain an integer number of bosons. Summariz- 
ing, the values of a and b are expected to be constrained 
by the requirements that 



q q ab 



are all integers. 



(2.39) 



insulating states with a and b as small as possible. Moti- 
vated by the application to the cuprates at hole density 
S = 1/8, let us look for solutions in which 6 is 4 or smaller. 
From Eq. (1.15) it is clear that a simple way to achieve 
this is to take tpg — unless ^ is a multiple of 4, so that 
Pmn is zero unless n is a multiple of 4. A generic solution 
with this property has p„in non-zero for all m, and hence 
a period a = 16. To obtain a smaller value of a we must 
make the pmn vanish for as large a set of m values as 
possible. However, with this particular choice of the cpg 
notice that p^n is proportional to Pm+4,n- So the best 
we can do is to have pmn non-zero only for m a multiple 
of 4. In this case we obtain a = 4, which results in a unit 
cell consistent with Eq. (2.39). 

We used the above strategy to identify a number of 
tpi configurations with highly symmetric, small unit cell 
density modulations. A few sample results are shown in 
Fig 8. All unit cells are clearly consistent with Eq. (2.39). 



E. Renormalization group analysis 

We now briefly address fluctuation effects across the 
mean field transitions found in Section II D. The anal- 
ysis follows previous work'^^"'^^ on scalar quantum elec- 
trodynamics in three dimensions who examined the the- 
ory So + Si for the case where the quartic interactions 
had a full \J{q) invariance. Here we will extend the ear- 
lier one loop results to the particular quartic couplings 
in Eq. (2.21). Unfortunately, we find runaway fiows to 
strong coupling for all values of q we have studied, and 
no stable fixed points which are accessible in a one loop 
analysis. 

For the renormalization group (RG) analysis it is useful 
to write the quartic terms in Si in the more general form 



Irnni 



(2.40) 



where 

Uim;ni = 5i^ra.n+ilm~t,n-&- (2-41) 

These couplings obey the symmetry relations 



(2.42) 



(as always, all index arithmetic is modulo q) and the 
trace identity^° 



Ulrn-li = US„ 



(2.43) 



To illustrate the origin of the constraints in Eq. (2.39), 
let us consider the g = 16 case, and attempt to construct 



with U = J2i 
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FIG. 8: A few tpf configurations for q — 16 with higher sym- 
metry and unit cells smaller than the maximal 16 x 16. All 
figures have ipe — unless ^ is a multiple of 4. (a) ipo = 1, 
</94 = 1, ips = 1, V'la = 1,4x4 unit cell. (6) Lpo — 1, ip4 — —1, 
ips — 1, (pi2 = 1, 16 X 4 unit cell, (c) (po — 1, ip4 = 1, 
(fis = —1, 'P12 ~ —1, 8x4 unit cell, (d) (po = 1, V'4 = 0, 
ips = 1, ipi2 = 0, 8 X 2 unit cell. 



The RG equations can be derived by a simple general- 
ization of the methods in Refs. 39,40. At one loop order 
the results are 



~dl 
dl 



Cq 



(e - 2ri)uim;nk - C 




2 '^£m;ij Uij-^nk 



7] = -aie^C. (2.44) 
The equations have been derived in 3 — e spatial dimen- 



sions and 7] is the anomalous dimension of the (pe field. 
The numerical constants C, ai_2 can be obtained both in 
the e expansion and in the fixed dimension RG performed 
directly in two spatial dimensions. For the e expansion 
we have = Stt^, ai — 3, a2 — 6. For the fixed di- 
mension expansion we have C^^ — Stt, ai — 8/3, a2 = 4. 
Using Eq. (2.41), the equation for the quartic couplings 
becomes 



d7„ 



dl 



(e - 2ri)j„ 




lmn2n-m.,n-l 



+ "yinjm+n-e,n + ll,m-n'~il.-n,m~n 



026 



'*Crf((5„o 







(2.45) 



We have verified that this flow equation for the couplings 
7mn preserves the symmetry constraints in Eq. (2.22) 
(for rectangular symmetry) and Eq. (2.24) (for square 
symmetry) . 

We searched for fixed points of Eqs. (2.44) and 
Eqs. (2.45) for small values of and found that all were 
unstable (as expected from Refs. 38,39). The analysis 
becomes rapidly more complicated with increasing val- 
ues of g, and we do not have general understanding of 
the fiow structure of these equations. 



F. Deviations from density p/q 

Our formalism has so far been explicitly designed for 
rational boson densities / obeying Eq. (1.11). In princi- 
ple, it would appear then we can model essentially any 
density simply by taking q large enough. However, this 
is not an attractive way to proceed for reasons discussed 
below. 

First, recall that there are q minima in the Brillouin 
zone, and so these minima become more closely spaced 
with increasing q. So, to consider these minima inde- 
pendent degrees of freedom (as is implicitly done in our 
continuum field theory), we can only work at momenta 
which are significantly smaller than the spacing between 
the minima. Hence the continuum theory of Section H C 
only applies beyond a large length scale which increases 
linearly with q. 

Second, the q minima are very shallow for large q. In 
other words, the bandwidth of the lowest lowest energy 
Hofstadter band containing these minima is very narrow. 
When integrating out the vortex modes away from the 
minima to arrive at the continuum field theory, one en- 
counters energy denominators determined by this band- 
width. In fact, this is likely the most severe restriction on 
the naive application of the q- vortex method. In particu- 
lar, one may argue that for p = 1 (and likely for any p by 
some refinement of this argument) the bandwidth of the 
lowest Hofstadter band is exponentially small, 0{e~^'^) 
with some 0(1) constant c, for large q.^^ This is a much 
smaller energy scale than the splitting between Landau 
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bands, which is also smaU but 0{l/q) for large q. A quan- 
titative physical consequence of the narrow bandwidth is 
a reduction of the vortices' tendency to condense, i.e. an 
enhancement of the domain of the boson superfluid state. 

As a consequence of these difficulties, it is preferable 
to take the continuum limit appropriate to the low en- 
ergy fluctuations of some nearby insulating commensu- 
rate state with a moderate value of q. The actual density 
of the boson system may not exactly equal p/g. However, 
it is easy to modify the theory for the system at density 
p/q (presented in Section II C), to allow for this density 
deviation. As the boson density is the average A^^ 'mag- 
netic' flux, we simply have to replace the action Sq in 
Eq. (2.19) by 

So = j d?rdT {y^ [\{d^ - iA^)^e\^ + s\^e?] 

+ ^ (e^^xd^Ax - ^ V) ) , (2-46) 

where Sf is the difference between p/q and the actual bo- 
son density, and a is the lattice spacing. All the higher 
order terms in the action remain as discussed in Sec- 
tion lie. 

The theory in Eq. (2.46) has a structure similar to 
the Ginzburg-Landau model for a 'superconductor' in a 
'magnetic' field, and has corresponding possibilities for 
its phases. It can behave like a type I 'superconduc- 
tor' and expel the flux df, while condensing the ipi: this 
clearly yields the commensurate Mott insulator with den- 
sity p/q. In the 'normal' phase of this 'superconductor' 
the ipi are uncondensed and flux df pierces the system: 
this is the superfluid with density p/g'-|-^/. This 'normal' 
phase can also have a condensate of vortex-anti-vortex 
composites with some of (pmn) non-zero, and this is a 
supersolid with density p/q + 6f, but with density wave 
order similar to that in the commensurate Mott insula- 
tor. Finally, the model Eq. (2.46) can also behave like a 
type II 'superconductor' and allow the 'flux' to penetrate 
the system in an Abrikosov lattice: this is an incommen- 
surate insulator (or floating Wigner solid) in which the 
particle density is p/q + Sf. Partial 'flux' penetration 
is also possible, and this yields Mott insulators at other 
densities. 



III. BOSON FRACTIONALIZATION 

The discussion of the previous sections generalizes the 
dual effective action of Refs. 15 (and its predecessors) 
for / = 1/2. In that work, it was pointed out that for 
/ = 1/2 the theory with 701 < is itself equivalent to 
a fractionalizcd representation in which the elementary 
boson is split into two "half" bosons, carrying equal (^ 
the elementary value) physical charge but opposite val- 
ues of an emergent U{1) gauge charge. These fractional 
particles are deconEned at the critical point in the sense 



that the gauge magnetic flux is conserved by the fixed 
point theory. In the language of the present paper, this 
conservation law is nothing but the conservation of the 
difference of the two (Co and Ci) vorticities. 

Two important features are critical to this conclusion. 
First, we used the equivalence of the dual 2— vortex ac- 
tion - in a particular regime (701 < 0) to a fraction- 
alizcd 1/2-boson representation. This mapping is non- 
trivial and indeed does not generalize to all /, as we 
will discuss below. Second, even when such a representa- 
tion does obtain, deconfinement at the QCP is contingent 
upon the fixed point theory possessing the appropriate 
emergent conservation law(s). This is true for / = 1/2 
roughly because the most relevant operator to violate the 
conservation law ((CoCi)^ + li-C-) is sufficiently high order 
that it can be persuasively argued^^ to be irrelevant. In 
general, the relevance or irrelevance of such operators is 
far from obvioTis, as the interacting g-vortex gauge the- 
ory in 2 + 1 dimensions is strongly coupled, so that the 
critical indices needed to answer this question are not 
calculable analytically. For this reason, we will have lit- 
tle to say in this paper about this second issue, which is 
best resolved by numerical (Monte Carlo) studies. 

In this section we will however discuss the first require- 
ment for a fractionalized boson interpretation of these 
Mott QCPs. Perhaps surprisingly, we will see that this 
can be satisfied in a theory with full square lattice sym- 
metry only for a rather special set of values of ratio- 
nal /. We do not at present have a general criteria 
to determine these values of /. By explicit example, 
we find that a fractionalized formulation is possible for 
/ = 1/2.1/4,1/8,1/9, and have proven (see below and 
Appendix E) that none exists for /= 1/3. We speculate 
that the sequences of / = p/q with q = n'^, 2r? and in- 
teger n all admit a fractionalized formulation. If the mi- 
croscopic symmetry is relaxed from square to rectangular 
(i.e. invariance under only twofold rather than fourfold 
rotations), however, the fractionalization requirements 
can be satisfied for any rational /. We show how these 
conclusions can be obtained both from an analysis of the 
dual g-vortex action of this paper, and from a comple- 
mentary direct analysis of the original boson problem. 
The agreement between these two completely distinct 
methods provides an important check on the results, as 
well as useful physical interpretation. The direct method 
is a generalization of the J7(l) gauge theory approach 
employed in Ref. 15. 

Although this analysis provides a formal means of un- 
derstanding fractionalization at such Mott QCPs, a more 
physical picture is also desirable. In Ref. 45, Levin 
and Senthil provided such a simple real-space picture of 
charge- 1/2 excitations for / = 1/2 coming from the Mott 
side of the transition. They showed that these particles 
can be understood as vortices points of connection of 
four elementary domain walls - in the density wave (va- 
lence bond solid) order parameter. For the specific case 
of / = 1/4, we discuss, in Section IIIC, a similar picture 
- as Z4X Z4 vortices - for the charge 1/4 excitations that 



16 



(may) become deconfined at the Mott transition. 



A. EVactionalization in the vortex picture 

Let us consider the general g-vortex action at fiUing 
/ p/?- By blind analogy with half filling (/ = 1/2), 
one might expect this to describe a QCP in which the 
elementary boson has been fractionalized into q bosons 
with charge 1 /q. If wc examine this proposition in a little 
detail, it becomes clear that this is not quite generally 
correct. 

For half- filling, -'^^ this conclusion was based on a dual- 
ity mapping of the 2- vortex theory back to a half-boson 
theory. The physical content of this mapping is that a 
'vortex' in a single vortex field ipi (this is a sometimes 
confusing double duality which unfortunately cannot be 
avoided in this particular method of analysis) is a point 
particle whose creation or annihilation operator becomes 
the fundamental field in a new representation. Here and 
in the following wc use single quotes to distinguish this 
'vortex' from a true physical vortex, which is created by a 
linear combination of ipi fields. As in any other Ginzburg- 
LandaTi-likc theory, such a 'vortex' will bind a quantized 
amount of 'flux' JcPr eijdiAj. The 'flux' however is 27r 
times the physical charge (boson number), so this 'vor- 
tex' is actually a particle carrying some amount of boson 
charge, which we will see is fractional under favorable 
circumstances. 

The above argument has one crucial un-stated assump- 
tion: that vortex-like topological defects in the fields 
are well-defined. This is true only if the magnitudes \(pi\ 
are constrained to be approximately constant, i.e. that 
the dominant terms in the action are minimized by this 
condition. Moreover, phases of these fields should be un- 
constrained. Formally, this requirement is necessary to 
pass from the soft-spin Ginzburg-Landau theory of the 
(fie fields to a hard-spin description in terms of phase- 
only variables. Duality transformations are formulated 
in a hard-spin description. 

Actually, we can relax this requirement slightly, as al- 
ready hinted by the mean-field analysis of the previous 
section. In particular, since the quadratic action Sq is 
invariant under a global U{q) unitary rotation of the 
(fie fields, it is sufficient that some unitarily transformed 
fields Q = Uu'ifie' (with V'^U = 1) can be regarded 
as having approximately constant magnitude, \Q\ = ^, a 
constant, but arbitrary phases. One may think of this as 
follows. Consider the evolution of the system from the 
superconducting to insulating state, by reducing s start- 
ing from s > 0. Due to fluctuations, the actual location 
of the critical point s = Sc &t which vortex condensation 
occurs is reduced to Sc < 0. Thus one may expect that 
locally, for Sc ^ s < 0, the Q flelds may develop some 
non-zero magnitude. This occurs because of a balance 
between the negative quadratic term — |s|^^|C^p and 
the quartic terms in £4. From the form of the negative 
quadratic term, this balance clearly favors states with 



a non-zero generalized "radius" yJ^TK^P in Ce space. 
The preferred directions in this 2g-dimensional (due to 
real and imaginary parts of each Q) space, however, are 
determined by the terms in £4. To obtain a phase-only 
description, one needs that the dominant terms in £4 
favor states with \Q\ =lp. 

This requirement places strong constraints on the form 
of £4. In particular, one must have that, for fixed 
^e IC^P — 9^^' minima of £4 satisfy = ^ for 
all £. Explicit examples with rotational invariance can 
be seen for g = 2,4 from Sec. IID1,IID3. For g = 2, 
this is true if 701 < 0. For g = 4 this is true in region 
E when A = 0. Thus if A is irrelevant in region E at the 
QGP, the requirement is satisfied. On the other hand, 
for (7 = 3, it is clearly not satisfied. For 701 = 0, there 
is a full 2g — 1 = 5-sphere of degenerate states, clearly 
of difiierent topology than the 3-torus of states of three 
angular variables. For 701 > 0, the degeneracy of min- 
ima of £4 is just a discrete multiple of the overall U{1) 
(singular angular degree of freedom) gauge symmetry. 

We do not have a general analysis of when the con- 
stant amplitude condition can be satisfied. One can, how- 
ever, place a necessary condition on the existence of any 
Lagrangian invariant under the symmetries of the prob- 
lem for which all minimum action configurations satisfy 
\(f \ =Tp and conversely all such configurations have min- 
imal action. In particular, for any such configuration, 
any new configuration equivalent to it by symmetry must 
have equal action. Thus if these are the unique configura- 
tions with minimal action, all symmetry operations must 
preserve this condition. That is, the space of states with 
\C,l\ ='(fi should be closed under all symmetry operations. 
Unfortunately, this is not trivial to check, since we still 
have a large freedom to choose the unitary transforma- 
tion Uui defining the C,e fields. 

For the simpler case in which the requirement of ro- 
tational invariance is removed, i.e. for a system with 
rectangular rather than square symmetry, it is straight- 
forward to see that this condition can be achieved for all 
/. This is easily seen from Eqs. (2.21) (2.22). In particu- 
lar, a quartic term of the form £4 = | \ is allowed 
for all q, and favors equal amplitude states. It is clearly 
invariant under all symmetry operations save rotations. 
Hence no unitary rotation is required, and we can simply 
take Q = V>t- 

The case with fourfold rotational symmetry remains 
an interesting open problem. While we have not yet 
obtained a general solution, we formulate the problem 
mathematically. In particular, the requirement that the 
constant magnitude configurations be closed under the 
action of the PSG is quite strong. In particular, an ar- 
bitrary such configuration, 1^^ = Ipe'^'^'- , must map to an- 
other such configuration. Since the transformations are 
linear, and the -di are arbitrary, this requires that the 
representation matrices of the PSG elements in this ba- 
sis must take the form G = AP, where, as we noted in 
Section I, A is a unitary diagonal matrix (i.e. a matrix 
whose only non-zero elements are complex numbers of 
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unit magnitude on the diagonal), and P is a permuta- 
tion matrix. Thus it is sufficient to determine whether 
the algebra of the generators of the PSG (unit x and y 
translations, and tt/2 rotations) can be faithfully repre- 
sented by such matrices. 

We believe that such permutative representations of 
the full (square symmetry) PSG exist only for specific q. 
As indicated above, we have constructed explicit exam- 
ples in detail in the text for g = 2,4. The ability to do 
so in these two cases is not obvious in the Landau gauge 
construction of the previous Sections. In Appendix C, 
we show that it can be clarified for 5 = 4 (for g = 2 it is 
rather straightforward) by choosing a different "symmet- 
ric" gauge which more closely respects the fourfold rota- 
tional symmetry of the lattice. The dual vortex theory 
constructed in this gauge immediately yields a permuta- 
tive representation of the PSG. Generally, a permutative 
representation exists only for specially chosen q. In Ap- 
pendix E, we prove that no such representation exists for 
q ~ 3, and also give explicit examples of such representa- 
tions for q = 8,9. The latter examples are obtained from 
a construction which can be generalized to g = n^,2n^ 
with any integer n, but which we cannot guarantee yields 
a solution. Nonetheless, we speculate that all members 
of this family support a permutative representation. 

For those cases for which the phase-only description 
can be achieved, we return to the question of charge frac- 
tionalization by considering the 'flux' (charge) of the dual 
'vortices'. Supposing the constant amplitude condition, 
for slowly varying •&£ the effective phase-only action can 
be written 



dP'rdT 



EPs 
, 2 



1 
2? 



(3.1) 

with p'^ = 2</?^ . Consider a fixed 'vortex' (independent of 
r) in one - say i^o - of the q phase fields, centered at the 
origin r = 0. One has the spatial gradient Vi^o = 4>lr, 
all other i?£ = for £ 7^ (here — (— 2/:a;)/r is the 
tangential unit vector) . Clearly, the action is minimized 
for tangential A = A(\). Far from the 'vortex' core, the 
Maxwell term (e^^A^jyA^)^ is negligible, so one need min- 
imize only the first term in Eq. 3.1. The corresponding 
Lagrange density at a distance r from the origin is thus 

4 = §[(lA-A)' + (a-l)^']. (3.2) 

Minimizing this over A, one finds A = l/{qr). Integrat- 
ing this to find the flux gives 



A - dr = — . 



(3.3) 



Since the physical charge is just this dual flux divided 
by 27r, the 'vortex' in t^o indeed carries fractional boson 
charge 1/q. 

This simple analysis can be confirmed by a straight- 
forward but more formal duality calculation on a lattice- 
regularized hard-spin §1 model. In the continuum, the 



theory of the q ip^ vortex fields coupled to the single non- 
compact U(l) gauge field A^ is described by the action 
<So -|- / d^xdrCi where <So is as in Eq. (2.19) and £4 has 
the structure (in an appropriate basis which realizes the 
permutative PSG) 



(3.4) 



where the co-efiicients Vmn are constrained by the per- 
mutative PSG: specifically, the matrix v must be invari- 
ant under all the permutative transformations P of the 
PSG: V = P~^vP. As we have discussed above, this 
theory has q — 1 globally conserved charges associated 
with the symmetries of phase rotations of the differences 
of the 'd( fields. The naive dual form of this continuum 
theory obtained"^^'^^ from a lattice-regularized hard-spin 
model has the following continuum representation in 
a theory of q bosons S^g (these are the 'vortices' in the 
vortex field (pe) each carrying boson charge 1/q, and q 
non-compact U(l) gauge fields A^ with the action 

fipa^pA^ j 



\e=o 



(3.5) 



As below Eq. (3.4), the matrices K and v are also in- 
variant under the permutative transformations P of the 
PSG. Note that the 'center-of-mass' gauge field Aj^ is 

'Higgsed' out into a massive mode by the coupling K. It 
is therefore convenient to make a new choice of variables 
in which this 'center-of-mass' mode has been explicitly 
set to zero 



A^' — A^' - A^~ 



(3.6) 



In these new variables, the J^e A^^ gauge field is decou- 
pled from all matter fields, and so is a redundant degree of 
freedom which drops out. The conserved gauge fluxes of 
the remaining (g— 1) gauge fields represent the (g— 1) 
conserved currents of the global symmetries of the ip^ 
theory noted above. The field theories discussed in this 
paragraph have a convenient representation in 'moose' or 
'quiver' diagrams"*^^, as illustrated in Fig 9. 

The following subsections will present a careful analy- 
sis of the reverse of the above duality, with special atten- 
tion being paid to lattice Berry phases and symmetries. 



Practionalization in the direct picture: U{lY ^ 
slave bosons 



Having seen that the dual g- vortex theory is under par- 
ticular circumstances itself dual to a theory of g charge 
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FIG. 9: Graphical representation of the field theories dis- 
cussed towards the end of Section III A in "moose" diagrams 
for q = 4. The circles represent non-compact U(l) gauge 
fields, while the directed lines represent complex scalar fields 
which are charged under the gauge field site at their 

beginning/end. The ifie vortex theory in (a) is specified in 
and above Eq. (3.4) with = Aqu — Ai^. Its dual theory 
of charge 1/q bosons (which are 'vortices' in ipi) in (b) is 
specified in Eqs. (3.5) and (3.6). In both diagrams the center- 
of-mass gauge field, Ao/j, + Ai^ or , decouples from all 

matter fields, and is therefore immaterial. If we visualize the 
Ao/Ai gauge fields on the north/south pole of a sphere, and 
the ring in (b) along the equator, then the geometric inter- 
pretation of the duality becomes evident. It is also amusing 
to note that the two diagrams are identical only for q — 2, 
and hence only this theory is self-dual, as noted by Motrunich 
and Vishwanath.*'* 



1/q bosons, it is natural to ask whether this fraction- 
ahzed representation can be directly obtained from the 
original boson problem, without recourse to duality. We 
show how this is done here, and further how the g- vortex 
theory can be recovered from this route by dualizing the 
fractional boson theory directly. 



TT-ii = n-i. In this subsubsection we will keep q arbitrary, 
with no a priori connection to the filling / = p/q. This 
connection is made in the following subsubsection. 

Eq. (3.9) expresses q — 1 independent constraints, 
which generate the same number of independent unphys- 
ical gauge rotations of the (1)^. That is, any rotation 

(t>ii ~^ 'Pil + Xii with Xii — leaves the physical bo- 
son creation/annihation operators invariant. It should 
be noted that the representation above has another set 
of discrete (Sq) gauge symmetries: the q slave rotor vari- 
ables can be permuted arbitrarily and independently on 
each site without changing fii or (pi. If the constraint in 
Eq. (3.9) is imposed exactly on every lattice site, how- 
ever, the Sq gauge symmetry does not lead to any addi- 
tional constraints, and it leaves the U{1)^~^ constrained 
states invariant (i.e. these states form a trivial identity 
representation of Sq). 

Such a microscopic slave-rotor rewriting of the problem 
has, unfortunately, no physical content. It becomes phys- 
ical only if, on long scales, a set of renormalized effective 
degrees of freedom "descended" from the slave rotors be- 
come physically meaningful quasiparticles, which locally 
violate these hard constraints. We would like, therefore, 
to consider a renormalized effective theory in which the 
fractional rotors have some local independence. In prin- 
ciple such a theory can be constructed by implementing 
the hard constraints in a path integral formulation, and 
studying fluctuations around a mean-field solution of this 
lattice field theory. For the universal phenomenological 
purposes of this paper, however, this is unnecessary, and 
we instead prefer to proceed on conceptual grounds. Im- 
portantly, it is essential that globally, in any finite sys- 
tem, all quantum numbers be physical, i.e. the total 
number of bosons must be integer. This goal is satisfied 
along with local fluctuations of the fractional rotors by 
introducing dynamical compact U{1) gauge fields: 



1. Construction of an effective gauge theory Hamiltonian 

Starting from a microscopic Hamiltonian like Eq. (2.2), 
one may introduce charge 1/q slave-rotor variables, 
(pUTfiie, with [(t)ii,hjef] = iSij5u>: 



9-1 

1 



q 



(3.7) 
(3.8) 



Clearly, the operator nu counts the number of flavor £ 
bosons, each of which has physical charge 1/q. This is 
a microscopically faithful representation of the original 
rotor states if we, e.g. impose the constraints 



flu - hii' = 



(3.9) 



for all £,£'. Thus, the physical state with fii = is 
represented by the state with all slave rotors identical 



(3.10) 



with {a,f3} e {x,y}. Here has integer eigenvalues, 
and Afp is a 27r-periodic angular variable. To keep the 
formalism as symmetric as possible, we have introduced q 
{£,£' = .. . q—1) such fields despite there being only q— 1 
independent constraints. The redundancy is eliminated 
by requiring 
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The vector overline indicates the usual two-dimensional 
vector {x,y) notation. With these fields, we replace 
Eq. (3.9) by a set of Gauss' law constraints 



A- 5* 



(3.12) 



where, as usual, 



and the A- symbol indicates 



a lattice divergence. If we do not allow the gauge elec- 
tric flux to exit the system boundaries (or there are no 
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boundaries), Gauss' theorem applied to Eq. (3.12) im- 
plies as desired that the total boson charge in the system 
(4 ^i^nu) is integral. The associated effective Hamilto- 

nian is then naively 

u 

+1 E l^'l' - ^E "'^"^"(^ ^ ^')- (3.13) 

it at 

In the second line above we have added quotes around the 

cosine term to indicate that care must actually be taken 
here to ensure that the unphysical center of mass of the 
gauge fields (^^ A^^ is decoupled from the other physical 
orthogonal linear combinations. This will bo achieved in 
the path integral formulation to which we will soon turn 
through the use of a Villain potential. 

We now come to a subtle and crucial point in the 
gauge theory formulation. Eqs. (3.12), (3.13) incorpo- 
rate on long scales the gauge constraints necessary to re- 
move the unphysical aspects of the [/(1)^~^ gauge redun- 
dancy. Further, the Gauss' law relations as formulated 
in Eq. (3.12) largely remove the Sq permutation gauge 
symmetry. What remains of the latter is an invariance 
under some group of global permutations, under which 



(3.14) 



where P{i) is a permutation on ... g — 1. These permu- 
tations P{£) are not arbitrary: only for particular choices 
can linear transformations of A^ and £^ be chosen along 
with Eqs. (3.14) to keep Eq. (3.13) invariant. This global 
permutation invariance in general depends upon the par- 
ticular value of q. At a minimum, it includes all cyclic 
permutations (generated by the elementary cyclic permu- 
tation P{e) =1 + 1, for which A'^ ^ and £^ £^+^ 
leaves Eqs. (3.13) invariant), but it is larger for special 
values of q. The permutation symmetry has been re- 
duced to a global one since the gauge fields couple dif- 
ferent lattice sites, so bosons can no longer be permuted 
locally. This permutation invariance, if the microscopic 
constraints, Eq. (3.9), are strictly enforced, has no conse- 
quence. In our effective theory, however, for which states 
violating these constraints are included, this global sym- 
metry has meaning. Since the original microscopic boson 
model has no corresponding symmetry, this will lead to 
unphysical aspects of the effective theory if it is not cor- 
rected somehow. 

An appealing and simple way to remove the permuta- 
tion invariance is to add a term to the Hamiltonian which 
removes this symmetry. We will construct such a term, 
guided by a few requirements. First, it should reduce to 
a trivial constant if the hard constraints in Eq. (3.9) arc 
enforced exactly. Second, it should not affect the aver- 
age density of the slave bosons. Third, for reasons that 
will become apparent in the next subsubscction, it should 
enlarge the naive unit cell of the fractional boson kinetic 



terms. Fourth, it should allow for some projective imple- 
mentation of all the symmetry operations of the lattice. 
The last requirement is the most non-trivial, and more- 
over seems somewhat at odds with the third one. As we 
will see, it precludes us from constructing an appropri- 
ate slave-boson theory at an arbitrary filling factor. Our 
tentative choice is to include a term of the form 
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(3.15) 



Here rrij^ is an integer field defined by introducing q sub- 
lattices (the choice of sublattices will be made later), with 



ruie 



for i e sublattice 
otherwise 



(3.16) 



This term clearly satisfies the first three requirements, 
but the implementation of symmetry operations is prob- 
lematic. In particular, as we have actually tried to do 
(requirement 3), Eq. (3.15) breaks the naive translational 
and rotational invariances of the original model. Indeed, 
a non-trivial implementation of these symmetries is possi- 
ble only for special values of q and/or lattice symmetries. 

A relatively simple case is that of rectangular symme- 
try. In this case, it is natural to choose the set of q sub- 
lattices by defining the i?*'' sublattice as the set of sites 
i for which = £(mod q) - just labelling sequentially 
columns by £ = 0,l,...g' — 1,0, — Then translations 
along y are trivially preserved, while an x translation is 
implemented by translating by one lattice spacing and 
simultaneously cyclically permuting £ ^ £ + 1. 

Unfortunately, this prescription fails to provide a 
means of implementing the 7r/2 rotation operation of the 
square symmetry group. In general, we have been unable 
to find any choice of sublattices which does so for arbi- 
trary q. In fact, the only two cases for which we have 
succeeded in finding a sublattice choice that fully imple- 
ments the square symmetry group are q = 2,4. For q = 2, 
it is achieved by taking the two checkerboard sublattices 
of the square lattice. Then both x and y translations 
and tt/2 rotations about a dual lattice site must be ac- 
companied by a (there is only one!) cyclic permutation, 
while a 7r/2 rotation about a direct lattice site requires no 
permutation. For g = 4, one defines the four sublattices 
according to 



i e sublattice < 



for {i,,iy) = (0,0)(mod2), 

1 for {i,,iy) = (l,0)(mod2), 

2 for (ix, iy) = (1, 1) (mod2), 

3 for {ia:,iy) = (0,1) (mod2). 



(3.17) 



These are constructed by choosing the sites of the 2x2 
square including the origin at its lower-left corner, la- 
beling them counter-clockwise 0,1,2,3, and translating 
this square vertically and horizontally to cover the plane. 
This can also be seen (less trivially) to preserve all sym- 
metry operations. For instance, a 7r/2 rotation about a 
dual lattice site should be combined with a cyclic permu- 
tation, and a unit translation in the x direction is effected 
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by 



Si 

£} 



fl-i+x,?. 



'-'t+x 



'^t+xi 
'^t+x^ 

'^l+X^ 



(3.18) 



and similar relations for (j)i£,Af. One can catalogue and 
verify the remaining symmetries for q = 4. 

Crucially, we see that it is not even possible in princi- 
ple to choose a set of sublattices that is preserved under 
rotations for most values of q, which is a minimal require- 
ment for successfully implementing rotational symmetry. 



2. Analysis of the effective gauge theory 

Having discussed the construction of the effective 
gauge theory Hq ^nl+n\ in Eqs. (3.13), (3.15), with 
the constraint of Eq. (3.12), we now turn to its analy- 
sis. Notably, up to this point, we have rather arbitrarily 
chosen to split the original boson into q parts, without 
any direct reference to the filling factor. From the dual 
analysis of Sec. Ill A, wo expect however that the phys- 
ically appropriate fraction corresponds to q = q. The 
connection is as follows. Due to the presence of the stag- 
gered potential /i^ , each of the fractional bosons moves in 
an environment of reduced translational symmetry, with 
a g-site enlarged Bravais lattice unit cell. At a filling 
/ = p/q, the average number of each fractional boson 
per enarged unit cell is then fq — pq/q. The smallest 
amount of fractionalization then consistent with an in- 
teger number of fractional bosons per unit cell is then 
q — q. Intuitively, with this condition, it is possible for 
the fractional bosons to form a featureless Mott insulator, 
apart from the effect of gauge fiuctuations. We therefore 
now focus on this case, q = q. The problem as described 
by Hq = Hl+n\ (Eqs. (3.13), (3.15)), is amenable to a 
simple direct analysis. 

Consider first the Mott phase. In Eq. (3.13), this oc- 
curs for largo u ^ t, for which nui will be approximately 
good quantum numbers with only only weak local fluc- 
tuations in the ground state. For a properly chosen 
range (fine-tuning is not necessary) of /is, one expects 
hu — pmu to be a qualitatively good approximation 
to the ground state. Having determined the state of 
the fractional bosons, one may then study just the pure 
U{\Y~^ gauge theory 



K^"cos"(Ax^'^)a, (3.19) 

at 



with the simplified Gauss' law constraint 



(3.20) 



where eu = p{mij - m^^^+i). 

In arriving at Eqs. (3.19), (3.20), we have, to a first ap- 
proximation, treated the Mott insulator as a dcconfined 
U{1Y~^ "spin liquid" with only gapped gauge-charged 
excitations. The particular liquid in question is charac- 
terized by the presence of q gapped boson excitations and 
the PSG of the previous subsubsection specifying the im- 
plementation of physical symmetries. As is well-known 
for q — 2, however, the dcconfined (or more precisely 
Coulomb) phase of a pure compact U{1) gauge theory 
in 2 -|- 1 dimensions is unstable to confinement via the 
proliferation of monopolc instantons, as first understood 
by Polyakov. This is true here as well, and the monopole 
unbinding will lead to the various possible density-wave 
orders associated with the Mott state. 

To see this explicitly, one can perform a simple duality 
transformation for the gauge theory. In its Hamiltonian 
version^^, this is simply the change of variables 



A X^^ 

27r 



2tt 



(3.21) 

(3.22) 



where £^ is an integer- valued classical c-number field sat- 
isfying 



(3.23) 



and Xi is an integer scalar field. Due to the redundant 

extra gauge field that has been introduced, we need to 
i 

enforce xt = £ = 0. The former can be im- 
plemented by including a large mass term, 'm^i^^Xl)'^ 
term in 7i, but for simplicity of presentation we will just 
keep this constraint implicit. One can verify that the new 
variables are canonical 



Xb 



(3.24) 



This reformulates the gauge theory in terms of gauge- 
invariant variables. The Hamiltonian becomes 



n 



C/(l) 



27rea^£:^|2 - ii:"cos"(27rB^) 



(3.25) 

To treat the theory analytically is simplest for large 
K ^ V, where fluctuations of are typically small. 
One may the naively expand the "cos" term. This ap- 
proximation, however, spoils the periodicity B^ <-!■ B'^ + l, 
or equivalently the 27rx integer constraint on xe- Hence, 
when making this expansion, one should include an ex- 
tra term in the Hamiltonian to favor integer values of xe- 
One obtains then 



K 



(7(1) 
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-Xcosxe 



(3.26) 
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with K = {2i:)^K,v = t;/(27r)^. It is useful now to de- 
— ( . 

compose E into transverse and longitudinal parts, 



2tt 



(3.27) 



Again, it is preferable to take 'Y^f^rji ~ 0. Provided one 
chooses 77f,7^ bounded at infinity, the latter is decoupled 
from xi. and contributes only an unimportant constant 
to 7i. Since xt has been promoted now to a continuous 
field, one can then make the shift xi ~^ Xi. ~^ Vt- O'cie 
obtains finally 



||A„Xa.P + y(i?^) 



at 

-\C0s{Xal + Vat) 



(3.28) 



This is just a set of q 2+1-dimensional lattice sine-Gordon 
theories (with the implicit constraint that the center of 
mass does not fluctuate, YlitXi = 0). To interpret the 
sine-Gordon terms, recall that and Xi are canonically 
conjugate. Hence, the operator e*-^'^ generates a unit shift 
of . It thus creates a single quantum of the gauge 
flux. Based on the corresponding space-time configura- 
tions, such operators are called "monopole" (instanton) 
operators. 

It is clear from this discussion that such monopole op- 
erators transform non-trivially under space group oper- 
ations. In particular, if for a particular space group ele- 
ment the site a g{a) and under the PSG £ G{t), 
one has 



al-Xa.1 



JXcit pi{Vg{a).G{e)-Va,e) 



(3.29) 



For large K v, the fluctuations of Xai are large 
on short scales, and one can integrate out these short 
scale fluctuations ("coarse grain") perturbatively in A. 
This gives a continuum theory in which the sine-Gordon 
terms are spatially averaged in a cumulant expansion. 
One expects on general grounds that all those and only 
those multi-monopole terms that combine to form scalars 
under the space group survive in the continuum theory: 



K 



(7(1) 
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^|A„X.P + f(sT 



+ V{{xi}) 



(3.30) 



where 



a,C 



Here the angular brackets indicate (cumulant) spatial av- 
eraging of the oscillating rjae factors over sites a of the 
dual lattice (it is sufficient in practice to average over 
the 2q sites in which ijat takes independent values). 

In general, this averaging leaves some non- vanishing 2- 
monopole and/or 4- monopole sine-Gordon terms A2,A4, 
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FIG. 10: Background electric field configuration for q = 
4. The other three background fields can be obtained by 
rotation. 
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FIG. 11: Background scalar fields -qi for q = 4. Rectangular 
unit cells are shown for i = 2, 3. The filled circle indicates 
the origin of the direct lattice. 



which depend upon the form of rjat- As Polyakov pointed 
out, although the fluctuations of Xat are large on small 
scales, they are always bounded in 2 + 1-dimensions, so 
that ultimately these non-vanishing sine-Gordon terms 
"pin" the xt flelds. This has the consequence of gap- 
ping out the "photons" of the putative Coulomb phase, 
and furthermore, giving a non-zero average to the single- 
monopole operators e'-^*. According to Eq. (3.29), such 
averages break space group symmetries, so the resulting 
state necessarily has some kind of density-wave order. 
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To make these manipulations concrete, wc give the 
example of / = 1/4, with the sublattice choice of 
Eq. (3.17). First, one must solve the Gauss' law con- 
straint, Eq. (3.23). Then, the fields r]ae are determined 
from Eqs. (3.27) by the Poisson equations, 

A'r/a^ = 27r (ca^gAaf ^) ^ , (3.31) 

and the requirement that they be bounded and satisfy 
J2e Vai = mod 27r. The results are shown in Figs. 10,11. 
For this set of r]ae, one can readily calculate the spatial 
averages m Eq. (3.30) up to O(A^). One finds the effective 
potential takes the form 

^({X^}) = - X] [•^^ cos(2x^ - ^) + A4 sinAxe 

e 

-A22 X cos 2{xe -Xe), (3-32) 
Kl' 

with A2 ~ A^, A4, A22 ~ A**. This indicates that for g = 4, 
individual monopoles occur doubled rather than quadru- 
pled as for q = 2, in the partition function. They are, 
however, created with the non-trivial phase factor above 
(7r/4 phase shift inside the cosine proportional to A2). 

Let us now turn to the superfluid phase, occurring for 
t S> M in Eq. (3.13). In this case, we expect the fraction- 
alized bosons to condense, i.e. neglecting gauge fluctu- 
ations, {e^'^') ^ 0. By the usual Anderson-Higgs mech- 
anism, the U{1)''~^ gauge fields will thereby acquire a 
gap. With all q fractional bosons condensed, of course 
the physical boson is also condensed, and one has a su- 
perfluid. 

What are the elementary excitations in this phase? 
One expects these to be; q flavors of vortices and anti- 
vortices, corresponding to phase windings in each of the 
q boson fields. Despite the presence of the gauge fields, 
these will have a logarithmic energy cost, since due to the 
constraint 'Y^nAi = 0(mod27r), the gauge fields cannot 
screen a single vortex. Thus from this (and indeed any 
other) universal point of view the superfluid is completely 
conventional. 

One can imagine, however, constructing "vortex exci- 
tons", i.e. a vortex in the boson field and an anti- 
vortex in l'^^ field. These excitations carry zero physical 
vorticity, and indeed can be screened by the gauge fields 
keeping the constraint satisfied. Most simply, a vortex in 
the and an anti- vortex in the (^-|- 1)*'* field is screened 
by a 27r (elementary) flux in the field . Thus each vor- 
tex exciton binds an elementary gauge flux. Conversely, 
an elementary flux of this type will necessarily bind a 
vortex exciton in the superfluid phase, as the "bare" 
gauge flux can only penetrate the superfluid in combi- 
nation with a vortex exciton (at finite energy) by the 
Meissner effect. 

Moreover, there is no microscopic conservation law for 
such gauge fluxes/vortex excitons, as they carry no net 
vorticity. Thus there will be some amplitude for tunnel- 
ing between states with different numbers of such objects. 



A space-time event at which such a tunneling occurs is 
nothing but a monopole of the type discussed in the pure 
gauge theory. We have already described how to calculate 
the amplitudes for such processes in the preceding discus- 
sion of the insulating problem. In the superfluid state, 
due to this binding of monopoles to vortex excitons, one 
may replace the notion of a flux-changing event by the 
creation of a vortex exciton. Equivalently, if one thinks 
of the action of such a term on a configuration with one 
vortex present, this may also be considered as a vortex 
flavor- changing event. 

Regardless of the language used, one may deduce from 
this discussion the form of the dual vortex action. In par- 
ticular, it should be written in terms of q vortex fields C,q, 
describing vortices in each of the q fractional bosons. By 
the usual considerations applicable in such a vortex de- 
scription, the vortex fields should be coupled to a single 
non-compact /7(1) gauge field, A^, whose 3-curl repre- 
sents the physical conserved boson current. As usual, 
in zero applied magnetic field, time-reversal symmetry 
implies symmetry between vortices and anti- vortices, so 
that the theory must be "particle-hole symmetric", i.e. 
relativistic in these variables. Provided the PSG con- 
tains enough symmetry between the different bose fields 
(as in all cases we have considered), the quadratic part 
of the vortex action is therefore fixed to have the form in 
Eq. (2.19). The fufl action 5 = 5o -t- Id^rdrViC^i) has 
an additional "potential" part. One expects 

nO)=^o(|OI) + ^A(0). (3.33) 

Here the first term is a vortex-conserving potential deriv- 
ing from terms not involving monopole events. It can be 
chosen as, e.g. Vb = "X^f lOl^i simply to favor equally 
all vortex flavors. The second term is deduced from the 
monopole terms already derived above by the rule 

^A(C£)=n{x4)le^x.^eC;C.,,' (3.34) 

the latter replacement encoding the coincidence of 
monopole with vortex flavor-changing events. Here c is 
a constant associated with amplitude of C,i. In the ex- 
ample of q = 4, one can readily see that the A2 term in 
Eq. (3.32) leads to a term in Vx gauge-equivalent to the 
A term in Eq. (2.36). 

Thus we have recovered the form of the vortex action 
studied in the remainder of the paper by this reason- 
ing. Of course, rather than the physical arguments used 
above, we could proceed rather more mechanically by 
performing duality transformations both for the [/(I) 
gauge flelds and for the fractional boson fields. This 
procedure recovers identical results, and is sketched in 
Appendix F. 

C. Approach to fractionalization from the 
density-wave state 

At half-filling, a very intuitive picture of fractionaliza- 
tion at the deconfined QCP could be developed based 
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upon topological defects in the valence bond ordered 
state. In this subsection, we consider the generalization 
of this notion to a potential deconfined QCP at / = 1/4. 
We will see that indeed there are intersections of par- 
ticular "elementary" domain walls in the charge ordered 
state that can be regarded as discrete vortices. Under not 
very restrictive conditions, such vortices trap a fractional 
charge in units of 1/4. Though we frame the discussion in 
this rather specific context, it is clear that the arguments 
in this subsection are rather general in nature. Indeed, 
we expect that it is a common (not universal but true in a 
sizeable fraction of cases) feature of charge ordered states 
breaking discrete lattice symmetries, that intersections of 
domain walls carry fractional charge. Of course, this does 
not imply in and of itself any exotic physics. Such dis- 
crete vortex configurations are certainly "confined" in the 
charge ordered state, having an energy linear in system 
size due to their trailing domain walls. 

Turning now back to the specific problem oi f — 1/4, 
we first need to understand the order parameter space of 
the density wave state. As discussed above, the appro- 
priate phase in question corresponds to that in region E 
of Fig. 5. From the figure, one can immediately see that 
the density wave has a 16 site unit cell with a period of 
4 lattice spacings in both the x and y directions. For 
reference, these have the following symmetry properties 
(see Sec. D): 

1. All 16 states are invariant under and Ty, i.e. 
the charge density pattern in the ground state has 
a period of 4 lattice constants in both x- and y- 
directions. 
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FIG. 12: Caricature of the ground state for the CDW phase 
at / = 1/4. Here the black circles are occupied sites. A unit 
cell is outlined in the lower left corner. Sites on the boundary 
are shared amongst 2 or 4 unit cells. The squares (open and 
filled) indicate the two inequivalent centers of 7r/2 rotations 
within each unit cell. 
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The unit cell contains two inequivalent centers of 
7r/2 rotation, separated by two lattice spacings in 
the X and y directions from one another. 



3. Refiection planes along the x and 
through the centers of rotation. 



y axes pass 



4. All 16 states are related to each other by successive 
Tx or Ty operations. 

All these symmetry properties, and the proper parti- 
cle filling, can be reproduced trivially by a wavefunc- 
tion taken as an exact boson number eigenstate on every 
site, of the form shown in Fig. 12. It is then clear that 
this phase may be interpreted as a charge density wave 
(CDW). Another caricature which is also useful to con- 
sider is a state which is a direct product of states defined 
on the 4x4 unit cells, with exactly 4 bosons per unit 
cell, but with charge fluctuations within each unit cell 
(see Fig. 13). More generally, the ground state of course 
contains further local charge fluctuations. Nevertheless, 
we expect that qualitative properties within the CDW 
phase can be determined by adiabatic continuity from 
the simpler caricatured states. 

Clearly, the order parameter space can be specified by 
giving the lattice coordinates of one of the two inversion 
centers, i.e. an ordered pair of integers (m, n) modulo 



FIG. 13: A less extreme cartoon of a unit cell of the 16- 
fold degenerate CDW ground state. Different shades of gray 
represent different boson densities. 



4. Let us now turn to topological defects. First consider 
elementary domain walls, which connect "neighboring" 
states in this order parameter space, i.e. where the two 
coordinates (m',n') and {m,n) differ by a unit transla- 
tion in x or y. Clearly, any one state can be separated 
from four others in this manner. A cartoon of such an ele- 
mentary domain wall, oriented for simplicity along one of 
the principle axes, is indicated in Fig. 14. An important 
question is the nature of the local excitations in the vicin- 
ity of the domain wall, assuming for simplicity that it is 
fiat (i.e. not in a rough phase) and pinned at the bound- 
aries, so that it has no translational degrees of freedom.. 
Theoretically, one may imagine that such a domain wall 
may be either insulating or conducting. In the simplest 
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FIG. 14: Cartoon of an elementary domain wall oriented along 
the X axis. Upon crossing from top to bottom, the state is 
shifted by one unit in the x direction. 



models, as in the cartoon of Fig. 14, the charge config- 
uration in the vicinity of the wall is fully rigid, and one 
expects it to be insulating. More precisely, the excitation 
spectrum remains gapped in the presence of the domain 
wall. Moreover, one can count the excess charge (above 
1/4- filling) per unit length associated with the domain 
wall. In the caricature of Fig. 14, the excess charge is 
clearly vanishing. Due to the gap in the excitation spec- 
trum, one expects this zero excess charge to be preserved 
as one perturbs around the caricatured state. 

Henceforth we will assume that the elementary domain 
walls are insulating and uncharged. Now consider the 
nature of vortices in the CDW order parameter. Some 
inspiration can be obtained from the point of view of 
the preceding sections. We have seen that charge ±l/q 
particles correspond to ±27r vortices in each of the (e 
fields. Thus one expects to find a total of 8 distinct such 
vortex-like excitations in the density wave order parame- 
ter. Indeed, by using the explicit transformations of the 
Q fields under translations in Eqs. (C5), one can straight- 
forwardly see that each such vortex corresponds to the 
intersection of four elementary domain walls. Across each 
such domain wall, the CDW pattern is translated to that 
of another of the 16 ground states. For instance, one 
finds that a positive vorticity vortex in Co corresponds to 
the sequence of operations: 



y X 



^ y^ X 



(3.35) 



More generally, upon moving clockwise around the center 
of the "vortex" in real space, one walks either clockwise 
or counter-clockwise in unit steps in the order parameter 
space. We denote these Z4 x Z4 vortices. A schematic 
illustration is shown in Fig. 15. 

Since both the domain walls extending away from the 
Z4 X Z4 vortex and the bulk CDW states are gapped Mott 
insulators, one expects this state to have a definite and 
quantized charge. Again, our cartoon wavefunction gives 
some insight. Both the bulk and domain wall wavefunc- 
tions can be written in a form with a definite 4 bosons 
per 4x4 unit cell. In the Z4 x Z4 state of Fig. 15, all 
sites are covered by such unit cells except one site in the 
center of the "vortex". This fact is inevitable by symme- 
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FIG. 15: A Z4 X Z4 "vortex" defect CDW phase. The site 
where 4 domain walls intersect does not belong to any unit 
cell and thus corresponds to an excitation of charge —1/4 (if 
empty as drawn). 



try: four domains related by unit translations that can 
be covered by non-overlapping unit cells are related by 
rotation around a site of the direct lattice. Thus one can 
write two such caricatured wavefunctions, having either 
zero or one boson on this central site. As the occupancy 
per site in the bulk is / = 1/4, such states have, relative 
to the ground state, a definite relative charge of —1/4 
(unoccupied central site) or -f3/4 (occupied central site). 
Because both the bulk and domain walls have been as- 
sumed insulating, we expect those quantized charges to 
be a robust feature of the Z4 x Z4 vortex, even for more 
realistic wavefunctions. This is because, although some 
charge may "leak" away from the central site, it cannot 
escape through the insulating environs, so the total ex- 
cess charge summed over the region around the vortex 
will not change. 

It is notable that one does not obtain in this way two 
symmetry-related states with charges Q = ±1/4. The 
difference from the case of / = 1/2 is that there is no 
particle- hole symmetric limit corresponding to / = 1/4. 
Indeed, there is no symmetry relating the charge Q = 
— 1/4 and q — -1-3/4 states, and generally these will have 
different energy. Of course, such a strongly localized form 
of the wavefunctions can be even approximately appro- 
priate only deep in the CDW state, where any particle- 
hole symmetry is very strongly broken. Thus there is no 
intrinsic contradiction with the expectations of the pre- 
ceding considerations of this section. Of course, deep in 
the CDW phase such a state has an energy linearly di- 
verging with system size, due to the four domain walls 
extending to the system boundaries. Its emergence as a 
deconfined "particle" in the spectrum near the critical 
point is not at all obvious from this perspective. 

Apparently if a deconfined QCP obtains in this case, 
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particlc-holc symmetry must be at least approximately 
restored in its vicinity, which is remarkable from this 
point of view. Its absence deep in the CDW phase, while 
not intrinsically troubling, makes an attempt to formu- 
late a theory of the QCP directly in terms of such micro- 
scopic vortex defects not very promising. We therefore 
leave it an open subject for future work. 



IV. IMPURITY PINNING 

This section will show how our approach of- 
fers a natural framework for describing recent STM 
observations^"^'^^ of modulations in the local density of 
states of the cuprate superconductors. As we noted in 
Section I and demonstrated in II, the bosonic theories 
presented so far apply also to paired electron systems, 
with the PSG determined in Eq. (1.6) by the density 
of Cooper pairs. We will interpret the observed density 
STM modulations as representing density wave order as- 
sociated with a proximate superfluid-to-insulator transi- 
tion. Indeed, as we have discussed in Section I, the PSG 
requires density wave order to be present near all such 
transitions in which the elementary vortices do not form 
composites. 

As we will demonstrate below, an important advan- 
tage of our focus on the tp^ vortex degrees of freedom, 
and their connection to density wave modulations, is 
that we are able to use the same theory to compute the 
impurity-induced pinning of density wave order in zero 
magnetic field, and also the vortex-induced density wave 
order in a finite magnetic field. At zero magnetic field, 
there are vacuum fluctuations of vortex-anti- vortex pairs 
which, when pinned by impurities, lead to modulations 
observable in STM. At finite magnetic field, there is a 
net excess of vortices over anti-vortices (say) but pre- 
cisely the same pinning potentials again lead to density 
wave order. This section will consider a simple toy model 
of the pinning process, with the aim of exposing the very 
general relationship between the zero field and finite field 
STM experiments. 

The perspective offered here differs considerably from 
previous analyses of the STM expcrimcnts^^^"''''' (although 
VBS order associated with vortices was suggested prior 
to the experiments in Ref. 56). In Refs. 46-48 the spin 
and charge order were treated as the primary fluctuat- 
ing degrees of freedom which were pinned by localized 
impurities. However, the connection of this pinning to 
the location of the vortices was indirect, and in a sense, 
had to be put in by hand. It was assumed that a large 
coupling W|Q| in Eq. (1.5) caused each vortex to induce a 
large pinning potential for the fluctuating density wave 
order. Our present approach dispenses with this ad hoc 
procedure, and shows how the connection between vor- 
ticity and density wave order emerges directly from the 
underlying quantum physics. 

We will work within the continuum field theory of low 
energy vortex fluctuations in the superfluid discussed in 



Section II C. As our purpose is mainly illustrative, we will 
only consider this theory within a leading Gaussian ap- 
proximation; our approach can, in principle, be extended 
to include the non-linear terms discussed in Section II C 
using the formalism presented below. However, the gauge 
potential has to be considered with some care, and 
cannot naively be dropped: we have to include an average 
value of Aj- which accounts both for the vortex "chemical 
potential" (proportional to the applied physical magnetic 
field) and for the interaction between vortices in the vor- 
tex lattice. We therefore consider the following action for 
the ipe vortices 



d'^rdT^^\{dr -iAr{r)) ipe\ + \daipe\ 



+ my Wei' 



(4.1) 



where the index a extends over the two spatial directions. 
This is clearly derived from the action in Eq. (2.19), writ- 
ten in a notation convenient for our purposes here: the 
energy gap towards creation of a vortex or anti-vortex 
is niyC^, and c is a velocity. Translational invariance is 
broken by adding a pinning potential: we assume a static 
impurity potential which couples linearly to the density 
operator: 



Sv 



j (frdr^Vrr, 



»WPmn(i',-r) 



(4.2) 



where pmn is specified as in Eq. (1.15). Since p^ni^i ''') = 
p_rn-n{r,T) wc must havc V^„(r) = V_™_„(r). We 
will now consider the properties of the action <So + Sy in 
zero and non-zero magnetic field in turn in the following 
subsections. 



A. Pinning in zero magnetic field 



In this case we can set At = 0, and perform a perturba- 
tive expansion in powers of Knn- A one-loop computation 
of the (pe susceptibility shows 



{Pmn{r)) = S{\Qimn\f j rfV no(r - r') F„„(r') , (4.3) 
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where the susceptibility is 



no(r) 



= -Q 
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So the pinned order decays exponentially with r over the 
length scale l/(2myc). 



B. Pinning in a non-zero magnetic field 

A non-zero magnetic field induces a finite density of 
vortices, and consequently we have to consider the theory 
in the presence of a non-zero vortex "chemical potential" . 
Here we wish to focus on the behavior of a single vortex 
localized near the origin, surrounded by an infinite vor- 
tex lattice. In this case, the quantity At has two distinct 
contributions: (i) a spatially independent vortex chem- 
ical potential, proportional to the applied magnetic 
field, and {ii) a r dependent potential created by the sur- 
rounding vortices of the vortex lattice. The latter poten- 
tial arises from the superflow around each vortex: in the 
present dual theory it is the logarithmic "Coulomb" inter- 
action energy associated with the gauge field created 
by the total vortex "charge" density. For the vortex un- 
der consideration, this "Coulomb" potential is assumed 
to have a minimum near r = 0, and we can safely per- 
form a quadratic expansion of the potential around its 
minimum. Collecting these contributions, we write 



1 



2„2 



(4.5) 



The second contribution above represents the interaction 
energy between the vortices, and the frequency ujy is de- 
termined by the superfluid stiffness of the bosons, and 
the spacing of the vortex lattice. At a more sophisticated 
level, we really have to consider the collective oscillations 
of the entire vortex lattice, and quantize the motion of 
all the vortices together^''; however, we will be satisfied 
here with the simple 'Einstein model' of the vortex lattice 
oscillation spectrum which is implicit in Eq. (4.5). 

A finite density of vortices requires that the vortex 
chemical potential be larger than the energy gap to cre- 
ate a vortex. Hence we have /i^, > rriyC^. In the 



analytic computations of this subsection, we will con- 
sider here only the additional density modulation cre- 
ated by the vortices, while neglecting the 'vacuum' quan- 
tum fluctuations of the vortex-anti- vortex pairs; provided 
\iAr — my(?\ <C m„c2, these 'vacuum' polarization con- 
tributions can be assumed to yield an additive contri- 
bution to the density modulation equal to that already 
computed in Section IV A. Technically, we can ignore 
the vortex-anti- vortex pairs simply by performing a low 
frequency expansion of Eq. (4.1) in the presence of a non- 
zero To this end we define 



1 



\f2AA._ 
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(4.6) 



and expand Eq. (4.1) to leading order in temporal deriva- 
tives and in \iAT — rUyC? \ . This yields the following famil- 
iar action for a bosonic vortex "particle" with q 'flavors' 



So 
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(4.7) 



We also have to insert Eq. (4.6) into Eq. (4.2) and so 
obtain the pinning potential acting on the vortices: this 
potential has a non-trivial structure in the 'flavor' space 
of the q vortices, and this will naturally be crucial in 
determining the induced density wave modulation. 

Just as in Section IV A, we proceed in a perturbative 
expansion in powers of Vmn- However, unlike the previ- 
ous zero magnetic field case, we now find that the pertur- 
bation theory is degenerate, and remarkably the leading 
order result turns out to be independent of the overall 
scale of the |ymn|- To zeroth order in the Vmn, notice 
that the vortex particle described by Eq. (4.7) has a q- 
fold degenerate ground state, with its spatial wavefunc- 
tion given by the ground state of the harmonic oscillator 



*^(r) = Ue ^ j exp 



myUjyr 



(4.6 



where Ue is, for now, an arbitrary unit vector in the inter- 
nal vortex space. The Ug vector is determined by diago- 
nalizing the pinning potential projected into this degen- 
erate ground state manifold. From Eqs. (4.2) and (1.15) 
we deduce that is the ground state eigenvector of the 
matrix Mcm, 



MimUm = £oUe, 



(4.9) 



where 
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Note that the overall scale of J^itm is immaterial at this 
order in perturbation theory because it has no influence 
on Lif. Finally, the density wave order associated with 
this vortex is determined by inserting Eqs. (4.6) and (4.8) 
into Eq. (1.15): 

(Pmn(r)) = — exp(-m„a;^r) 

X a;™"/2^W;ZY^+„a;^™. (4.11) 

e 

As promised, this result has no prefactor of Vmn '■ the de- 
pendence on the pinning potential is entirely in Eqs. (4.9) 
and (4.10) where the relative values of Vmn determine the 
orientation of the unit vector Ue in the complex q dimen- 
sional vortex space. 

We also note here that the above calculation can be ex- 
tended to allow for small deviations from the commensu- 
rate density of p/q. We do this by extending Eq. (4.1) as 
in Eq. (2.46), and include an average dual 'magnetic' field 
of strength 27r6f /a^ that acts on the vortices. However, 
it cannot be assumed that the dual 'magnetic' field is spa- 
tially uniform, and its space dependence must be com- 
puted self-consistently. We expect that the dual 'mag- 
netic' flux will be partially expelled from the vicinity of 
the vortices (where the dual 'matter' resides). In the 
direct boson language, this means that the density of 
bosons is not spatially uniform in the presence of the 
vortex lattice, and that the boson density is closer to p/q 
than average at the positions of the vortex cores. This 
effect has been discussed in Ref. 58 for the case of Mott 
insulators with integer filling. Let us assume that the ac- 
tual boson density deviation from commensurability at 
the vortex lattice positions is 5f, where |^/| < \5f\. We 
now need to determine the vortex wavevfunction by solv- 
ing the Schrodinger equation associated with Eq. (4.7) 
and a 'magnetic' field of strength 6f. This is equivalent 
to the problem of a particle in a parabolic potential and a 
magnetic field; the wavefunction of the vortex retains the 
Gaussian form in Eq. (4.8), and all subsequent expres- 
sions remain as in the paragraphs above apart from the 

shift in the frequency ojy ^ i'^l + (27r(5//(2TO„a^))^ j 



C. Discussion 

The above computations of pinned density wave order 
in zero and non-zero magnetic field have a number of un- 
usual features which deserve further comment. It would 
be useful to have tests of these features in future STM 
experiments. 

The pinned order in zero field in Eq. (4.4) has a conven- 
tional form similar to that obtained in other approaches. 
The induced order decays exponentially from the pin- 
ning site, on a scale 1 / (2m„c) related to the distance to 
the supcrfluid-to-insulator transition. This decay length 
equals the intrinsic correlation length for density wave 



order fluctuations in the superfluid, as might naively be 
expected. Also, for weak pinning, the amplitude of the 
pinned order is proportional to the strength of the pin- 
ning. 

In contrast, the finite magnetic field result in Eq. (4.11) 
has some remarkable features. Quite generally, for weak 
pinning, the envelope of the density wave order has the 
functional form of a Gaussian in the small r region where 
the harmonic approximation holds. The width of this 
Gaussian is not determined by the correlation length of 
the density wave order, but by the scale over which the 
vortex undergoes quantum zero point motion. The latter 
distance is determined by the inverse square root of the 
mass of the vortex, m„, and the vortex lattice oscilla- 
tion frequency, LUy. It should be possible to compute the 
magnetic field dependence of LOy , and thus quantitatively 
compare Eq. (4.11) with STM experiments in a varying 
magnetic field. Finally, the prefactor of Eq. (4.11) is in- 
dependent of the strength of the pinning field, but has an 
interesting dependence on lu-u, which places further con- 
straints on possible experimental comparisons; the inde- 
pendence on the pinning field strength can readily tested 
in numerical studies such as those in Refs. 36, 59, and 
60. 



V. CONCLUSIONS 

This paper has presented a general framework for de- 
scribing the superfluid-to-insulator transition on regular 
lattices. We examined simple model boson systems on 
the square lattice, but we expect our results to be quite 
general, and applicable also to paired electron systems 
exhibiting transitions between superconducting and Mott 
insulating states. The latter connection is established in 
more detail in a companion paper^'' II. 

There were three important themes underlying our ap- 
proach. The first is a feature common to several recent 
studies of quantum phase transitions with multiple order 
parametcrs^'''^''''^^'^''': the quantum mechanics of the de- 
fects of the order parameter in one phase induce correla- 
tions of the second phase. In our case, we focused on the 
vortex defects of the supcrfiuid state. These were found 
to have remarkable symmetry properties which linked 
them intricately to density wave order. The condensation 
of such vortices, apart from leading to the destruction of 
supcrfiuid order, also led to the simultaneous appearance 
of density wave order in the Mott insulator. This consti- 
tutes an important mechanism for the breakdown of the 
Landau- Ginzburg- Wilson theory of the phase transition. 

The second theme was the crucial role played by pro- 
jective representations of the space group of the lattice 
(the PSGs). PSGs have been highlighted in the recent 
work of Wen^'^ in his study of ground states insulating 
spin systems (equivalently, of bosons at half-filling). We 
argued that the vortices of the superfluid in the vicinity 
of a commensurate Mott insulator state transform under 
a PSG defined by Eqs. (1.8) and (1.10). The key defining 
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parameter of this PSG is the unimodular complex num- 
ber a;, and for a Mott insulator with average density of 
p/q bosons per lattice site a; was defined in Eqs. (1.9) and 
(1.11). Note that the superfluid (or a supersolid phase, if 
present) need not be exactly at this density, and our for- 
malism allows small variations in density away from that 
of the commensurate Mott insulator (see Section II F). 
The existence of such a PSG was shown to imply that 
there are q degenerate species of vortices in the super- 
fluid. Moreover, we argued that any linear superposition 
of such vortices was necessarily associated with density 
wave order at wavelengths which are integer divisors of 
q lattice spacing: consequently, as noted above, vortex 
condensation produces a Mott insulator with long-range 
density wave order. 

Section IV presented an analysis of the consequences 
of impurities in the superfluid. Any impurity breaks the 
lattice symmetry, and hence it also lifts the degener- 
acy between the q vortex species. It then follows that 
any localized vortex has an associated halo of density 
wave order, and we offered this as the fundamental ex- 
planation of the STM observations of Hoffman et al}^ on 
Bi2Sr2CaCu208+5. 

The third theme of our paper was one of particle num- 
ber fractionalization. This arose in an attempt to in- 
terpret the dual vortex theory of the LGW-forbidden 
superfluid-insulator transition in the language of the di- 
rect lattice bosons. Following the 'deconfined criticality' 
proposal of Senthil et al}^, we hypothesized a theory in 
which the boson fractionalized into q constituent parti- 
cles each with particle number 1/q. By performing a du- 
ality analysis upon such a fractionalized boson theory, we 
were able to reproduce the dual theory of q vortices, but 
only for special values of q. The selected values of q were 
those for which it was possible to obtain a permutative 
representation of the PSG. 

Moving beyond simple boson systems, our analysis 
makes clear that the key ingredient needed for other su- 
perfluids is the value of the parameter uj controlling the 
PSG of its vortices. With the knowledge of iv, the gen- 
eral structure of the field theory in Section II C is fully 
determined, and subsequent results follow. We invite the 
reader to now proceed to II, where we consider electron 
systems with short range pairing, and show that w is 
determined by the density of Cooper pairs. 
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APPENDIX A: VORTICITY MODULATIONS 

We have shown in the body of the paper that the vor- 
tex fields ipi are naturally connected to modulations in 
the generalized density (which could be any observable 
invariant under time reversal and spin rotations); explic- 
itly, the two were connected by the formula in Eq. (1.15). 
In this appendix we will examine the issue of the modu- 
lations in the vorticity. 

Following Ivanov et al.^^, we define the vorticity as 
a observable associated with each plaquette of the di- 
rect lattice (or each site of the dual lattice), equal to 
the sum of the particle (boson) currents flowing anti- 
clockwise around the plaquette. The vorticity is odd un- 
der time-reversal, and because all the states considered 
in the body of the paper are time-reversal invariant, its 
expectation value vanishes. However, when a magnetic 
field is applied, as in Section IV B, then the vorticity 
can have quite a rich spatial structure. In the superfluid 
phase, we know that the magnetic field induces a lattice 
of vortices. We also argued in Section IV B that there is 
a density modulation superimposed on this vortex lattice 
described by Eq. (4.11). Here we shall show that there 
is also a corresponding pattern of vorticity modulations 
at the wavevectors Qmn- Thus the appearance of "stag- 
gered vorticity" is quite a general feature of superfluids 
near a Mott transition, and not a special characteristic 
of certain spin liquids as implied by Ivanov et al.^^. 

We can obtain an expression for the vorticity cither by 
explicit derivation or by symmetry considerations alone. 
In the first approach, we notice that the vorticity is con- 
jugate to a magnetic field applied on the plaquette of 
a dual lattice, and the latter couples to the operator 
ip*dr'>pa — drtpafpa- In the sccoud approach, we need the 
action of the time-reversal operation, T, on the varies 
fields in the paper; a simple analysis shows that under T 
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From these mappings we can construct the simplest ob- 
servable which is odd under T, /^"^', Iy^'^\ and trans- 
forms like Eq. (1.16) imdcr translations and lattice rota- 
tions. We define Vmn to be the Fourier component of the 
vorticity at the wavevector Qmn defined in Eq. (1.14), 
and obtain (compare Eq. (1.15)) 

9-1 



fe—g:^-^'Pe+uj. (A2) 



29 



Here Sv{Q) is an unknown structure factor which de- 
pends smoothly on Q. 

Wc can now easily extend the discussion in Sec- 
tion IV B to obtain the mean vorticity modulation in the 
vortex lattice. The analog of the result in Eq. (4.11) is 



/T/ / x\ Sv {\Qmn\) mv(jJv / 2\ 



(A3) 



To discuss the translation and rotation operations in 
momentum space it is convenient to introduce the recip- 
rocal lattice to our triangular lattice. The reciprocal lat- 
tice is also a triangular lattice which, with respect to the 
original (dual) triangular lattice, is rotated by 30 degrees 
and spanned by the basis vectors 



bi = -I- -^e^ 



bo = 



x/3 ^• 



(B7) 



APPENDIX B: OTHER LATTICES 

We illustrate the extension of the PSG to other lattices 
by considering the case of bosons on a honeycomb lattice. 
We consider boson density / per site, and as in the case of 
a square lattice we write / = p/q, with p and q relatively 
prime integers and define to = e^"^^^ . The vortices reside 
on the dual triangular lattice, and so we have to consider 
the PSG of motion on a triangular lattice in the presence 
of a 'magnetic' field. 

The triangular lattice is spanned by the two basis vec- 
tors 



ai = 



a2 = -^e, 
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Here e^; and By are orthogonal unit vectors. Any position 
vector of the lattice can be written as a = aiai + a2a2 
with ai and a-i integers. It is also convenient to define 



ad = ai -I- a2 



(B2) 



Here we will use the Landau gauge A^t = Aai = and 
Aa2 = 2fai, = /(2ai + 1) . (B3) 

For the translation and rotation operators we obtain 
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As in the case of a square lattice, the translation opera- 
tors do not commute, and in addition to T1T2 = oj^T2Ti 
we have 



T1T2 = uj Td 



(B5) 



As expected, the following relations of the triangular lat- 
tice are satisfied. 
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We can now write any wave- vector k as k = fcibi -|-A:2b2. 
Using Bi ■ hj — 6ij we have k • a = fciai -|- ^202 and obtain 
for the translation and rotation operations in momentum 
space 
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where i' is an integer. Note that we find a different 
transformation form for the rotations depending upon 
the value of q (mod 4). There is also a compact expres- 
sion for the rotation operator valid for general q, but this 
obscures the underlying symmetries by containing terms 
which cancel each other in the case of even q. 

As in the case of the square lattice in Section II B, 
let us now consider the important vortex fluctuations at 
momenta where the spectrum has minima. If |A) is a 
state at a minimum of the spectrum with T2IA) = e^*'^^ ^ 
then acting successively with Ti on this state gives us 
degenerate states Tf\A) with T2 eigenvalue e"''^^^;-^^ 
This procedure actually does not completely span the 
degenerate manifold for q = Av -\-2, and we will limit our 
considerations for now to the other cases. 



1. q = 2v + \ 

For odd q the eigenvalues e~^^'^uj~'^^ are all different, 
and so the degeneracy of the vortex fields has to be at 
least q. For the Hofstadter problem on a triangular lattice 
we find that there are exactly q minima located at the 
wavevectors (0, Anip/q) with £ = . . .q—1. Analogously 
to the case of the square lattice in Section II B, wc can 
now proceed and find for the transformation laws of the 
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vortex fields ^pe, at these minima 



R. 



Ti 
Td 

dual 



(fe (feu' 



-2e 



-(2£+l) 



■7r/3 '■ 



p-i7r(g-l)/12 ^ 
m=0 



^/9 



(l-1m)l 



(B9) 



It can be verified that these transformation laws obey 
Eqs. (B5) and (B6). 



4z^ 



If (7 is a multiple of 4 we only have g/2 different vortex 

minima which arc located at the wavevectors (0, AiT(.p/q) 
with £ = 0...g/2 — 1. The transformations among these 
fields are 




odual 



These also obey Eqs. (B5) and (B6). 
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FIG. 16: Pictorial representation of the "symmetric gauge" . 
Arrows represent a factor e~"^^^. Magnetic unit cell labelling 
convention is also shown. 



eigenvectors are given by: 

v"" = (0, 1/2, 1/v^, 1/2) , 
= (-z/2,l/Al/2,0), 
= (l/V2,i/2, 0,1/2), 
= (1/2,0, 1/2, 1/x^), (C2) 

As before, concentrating on the long-wavelength vortex 
dynamics near the transition, the vortex fields can be 
written as linear combinations of the 4 low-energy vortex 
modes: 



APPENDIX C: VORTEX BAND STRUCTURE 
AND PSG IN THE "SYMMETRIC GAUGE" 

In this appendix we show how to directly obtain a per- 
mutative representation of the PSG for / = 1/4 by a 
judicious choice of gauge for the dual magnetic flux. We 
will loosely call this the "symmetric gauge". It is most 
conveniently expressed pictorially and is shown in Fig. 16, 
along with the corresponding magnetic unit cell labelling 
convention. A unique property of this gauge choice, that 
will play a crucial role in the discussion below, is that the 
Bravais nei of the magnetic lattice is invariant under 7r/2 
rotations. It is clear that such a gauge choice is not pos- 
sible in general, but is at least for filling factor / = p/q, 
satisiying the condition q = 2,n^,n = 2, 3, 4,... (for the 
square lattice). We do not, however, explore this beyond 
<? = 4. 

Directly diagonalizing the lattice vortex action in 
Eq. (2.8) we find 4 low energy vortex modes at wavevec- 
tors: 

ko = (0,-^2), 
ki = (^,-^2), 
k2 = (^,7r/2), 

k3 = (0,7r/2), (CI) 

which are defined in the first Brillouin zone of the mag- 
netic lattice and we have taken the magnetic lattice con- 
stant to be equal to one. The corresponding normalized 



(C3) 



t=0 



where R labels magnetic unit cells, a labels sites within 
each magnetic unit cell and 



(C4) 



Using Eqs. (2.8) and (C3) we can work out the trans- 
formation properties of the low-energy vortex modes Lp(^ 
under the symmetry operations of the microscopic Hamil- 
tonian (2.2). These transformations turn out to have the 
following form: 

Tx '■ Co ^ Ca, Ci C2, C2 — Cii Ca Coj 

Ty-- Co ^ Ci> Ci Co, C2 «C3, Cs K2, 
'■ Co Ci, Ci K2, C2 Cs, Cs Co, 

-^^2' • ^0 ~^ ^0' ^1 ~^ C2 K2, Cs Cl, 

If ■■ Co - icr, Cl - <o. C2 - C3, Cs - C2, 

y • Co ^ Cs , Cl ^ C2 , C2 ^ Cl , Cs ^ Co, 

• Co ^ -Co, Cl -Cl, C2 ^ C2, Cs ^ C3, 

/dual . /- , /- /-* /- /- 

y • Co Co > Cl -Cl , C2 C2 , Cs — * Cs • 

(C5) 

Here denotes a 7r/2-rotation about a direct lattice 

site, is a 7r/2-rotation about a dual lattice site, /^'^ 
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are reflections with respect to x,y axes passing through a 
direct lattice site, while I^"^^ are reflections with respect 
to 3;, y axes passing through a dual lattice site. Clearly, 
the above rotational and translational transformations 
are all of the form of a unitary diagonal matrix multi- 
plying a permutation matrix. Hence by composing each 
with a suitable global (gauge) phase rotation, they realize 
a permutative representation of the PSG. 

It is straightforward to see that the most general quar- 
tic potential invariant under these transformations has 
the same form as Eq. (2.36). Thus one directly obtains 
in this manner the end result obtained in the text by first 
using Landau gauge and then unitarily transforming to 
(i variables. 



APPENDIX D: SOME GROUND STATES AND 
THEIR PROPERTIES FOR q = 4 

In this appendix we present a few more details of the 
structure of the saddle points of effective potential in 
Eq. (3.34) for the equal amplitude state identified as E in 
Fig 5. Writing Q oc e^^^, the relevant term in the energy 
density (Euclidean Lagrange density) is (assuming first 
that A < 0): 

£4 = -|A|[cos2(i?o-'!?i)-cos2(i?i-i?2) 

+ cos2(i?2-i?3) + cos2(i?3-i?o)]. (Dl) 

In these variables it is not immediately obvious what the 
ground state is. We change variables to 



^0 = ^0 

1?3 



TT 

4' 



(D2) 



(D3) 
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^3 



and 



^2-^3 



■ 0,7r, 
: 0,7r, 

TT StT 

4'T' 
4'T' 

TT StT 

4'T' 



(D4) 



Thus the ground state has a 16-fold degeneracy. To see 
what do these states look like in real space it is sufficient 
to understand how they transform under the symmetry 
operations of the microscopic Hamiltonian Eq. (2.2). Us- 
ing Eq. (C5) one finds that the ground state manifold 
has the following properties. 

1. All 16 states arc invariant imder and Ty, i.e. 
the charge density pattern in the ground state has 
a period of 4 lattice constants in both x- and y- 
directions. 

2. There are 2 states that are invariant under rota- 
tions by 7r/2 about a given direct lattice site. These 
states are, in the original variables: 



and 



^1 - '4, ^1 



-?2 = ^,^2 



^3 = V (D7) 



In the new variables the mean-field energy density be- 
comes: 

£4 = -|A| cos2(^a?o-^?i-^) +cos2(i?i-i?2-|) and 
+ cos2 {§2 - ^3 - I) + cos2 (t?3 - i?o - I) 

It is now clear that the minimum is achieved when 

^a-h = 0,7r, 



These 2 states are related to each other by {T^Ty)^ 
transformation. More simply, this implies that each 
state has two inequivalent centers of rotation within 
each unit cell. 

3. The states are also invariant under and ly trans- 
formations through the centers of rotation. 

4. All 16 states are related to each other by successive 

or Ty operations. 

It is then clear that the ground state is a charge density 
wave (CDW), of the general form shown in Fig. 13. 

Now consider the case A > 0. In this case there are 16 
degenerate ground states as well. They are given by: 

~ - TT StT 

= 2'Y' 

- - TT StT 

= 2'T' 



1?0 - 



TT StT 

2'T' 



(D8) 



- - TT 37r 
^o-i?i = -4,^, 

- ~ TT StT 



^9, - 



(D9) 



(D5) 



TT 37r 

~4'T' 

These states have exactly the same symmetry proper- 
ties as the ones at A < and thus are physically the same 
states. The two ground states that are invariant under 

R2L in this case are: 
2 

^?0-^l = y,1?l-1?2 = y,1?2-1?3 = y, (DIO) 

and 

1?0-^l = 1?1 -1?2 = ~, 1?2-1?3 = (Dll) 
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APPENDIX E: PERMUTATIVE 
REPRESENTATIONS OF THE PSG 

As discussed in Sec. Ill A, the problem of determining 
whether a fractionaUzcd boson interpretation of the q- 
vortex theory at fiUing f = p/q reduces to a problem 
in representation theory. In particular, we must find q- 
dimensional matrices x, y, r representing the generators 
T^,Ty, R^'jf (for this append ix we neglect the inversion 

operations) 

X — XX ^ 

y = yy, 



r = rr. 



(El) 



where x,y,f are diagonal unitary matrices (i.e. diagonal 
matrices with all entries of absolute value one), and x,y, z 
are permutation matrices. We seek matrices of this form 
satisfying 



xy 



yr 

„4 



ujyx, 

y-\ 

x, 

1, 



(E2) 



with w*^ = 1 if and only if fc = O(mod q). The form of 
Eq. (El) requires that the permutation (hat) parts of the 
generators satisfy very similar relations to Eq. (E2), to 
wit: 



we have been able to prove there is no such representation 
for (7 = 3. This proof follows. First, note that f can con- 
tain either one 2-cycle or be trivial. Suppose first f = 1. 
Then from the third of Eqs. (E3), we see that x = y = p, 
some permutation which must be non-trivial according 
to the above consideration. But since f = 1, the sec- 
ond of Eqs. (E3) implies = 1. Thus, without loss of 
generality, we can take p = (12) (we use cycle notation, 
so this is just the permutation of the first two objects). 
Inserting this into xy = ujyx, one finds that, since the 
third element is not permuted, 533^33 = ^533^33, which 
is impossible. Thus our assumption that f = 1 must be 
false. 

Hence if such a representation exists, we must take f 
to be a 2-cycle, and without loss of generality, we may 
take r = (12). Thus P = 1, so from Eq. (E5), x"^ = 1 and 

= 1. Now at least one of x,y is non-trivial, so let us 
assume, again without loss of generality, that x ^ 1. It 
must be a 2-cycle since its square is unity. There are only 
two inequivalent possibilities. If a; = (12), then x = f, 
which implies y = (12) as well. We have already seen 
that there is no solution with x = y = (12), so this is 
inconsistent. The other inequivalent choice is i = (23). 
Then we can compute y = ■rxf~^ = (13). But it is easy 
to see that this violates the first of Eqs. (E3). 

Thus all assumptions for (7 = 3 have been proven in- 
consistent save one: that no such representation exists. 

We have found some explicit examples for = 8, 9. 

For q = S, with the notation 7 = e™^^^ (7^ = w), we 
find 



xy 



yr 



yx, 

X, 

1. 



(E3) 



This leads to a few general observations. First, recall 
that any permutation can be written (decomposed) as a 
product of disjoint cycles. Since = 1, f can contain 
only cycles of length 1,2, and 4. Hence, contains only 
1-cycles and 2-cycles, and an even number of the latter. 
Another interesting identity is 



2 2 
r yr 



— 2 2 

r yr 



Hence 



22 1 

r yr y = l, 



1 1 
r xr X 



1, 



(E4) 



(E5) 



the latter equation being obtained similarly. Thus 'Px 
can contain only 1-cycles and 2-cycles as well. 

Finally, to obtain a non-vanishing commutator of x 
and y, clearly at least one oix,y must not be the identity 
permutation. 

Unfortunately we have so far been unable to solve this 

representation problem generally. The explicit examples 
for (7 = 2,4 have been discussed in the text. In addition. 
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For q = with the notation uj = e^'^/^ (d)^ = w), we 
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find 
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These were obtained as follows. First, we obtained 
i, y, f from a simple physical construction. In particular, 
we constructed a covering of the sites of the dual square 
lattice by q sublattices. For g = (n = 3 for g' = 9 
above), these are chosen by breaking the lattice up into 
a square grid of n x n blocks, labeling the sites in each 
block identically from £ = \ . . .q. Taking the site of the 
direct lattice at the center of one of these blocks as the 
origin, one can see that rotations and translations all act 
as permutations of the q labels. These pcrmiitations arc 
taken as x, y, f and by construction satisfy Eqs. (E3). For 
q = 2r? (n = 2 for = 8 above), we break the sites of the 
dual lattice again up into a square grid of n x n blocks, 
and further break up these blocks into staggered A and 
B sublattices. We label the sites of each block of the A 
sublattice identically from \ . . .v? and the sites of each 
block of the B sublattice identically from n'^ + 1 . . . 2r? . 
Now taking the origin as a point of the direct lattice at 
the corner of two A and two B blocks, we once again 
obtain from physical rotations a set of permutations on 
q variables, which solve Eqs. (E3). 

We then take x,y,r as unknowns, and look for solu- 
tions of Eqs. (E2). These reduce to many redundant 
equations for 3q unknowns (the phases of each diagonal 
entry ofx,y,f). We find for g = 8,9 many solutions, of 
which the above examples are particular ones. This pro- 
cedure can clearly be generalized to arbitrary n as defined 
above. We do not, however, have a proof at present that 



solutions for the diagonal unknowns exist in general. The 
overabundance of solutions for g = 8, 9 suggests, however, 
that they do exist - hence the speculation in the text. 



APPENDIX F: DUALITY FOR g-SLAVE BOSON 
GAUGE THEORY 



In this appendix, we directly dualize the hamiltonian 
Hq = Hl+n\ in Eqs. (3.13), (3.15). The first step is to 
shift the slave rotor number operators to new variables 
with an average zero density: 



nu flu +pmu. 



(Fl) 



Note that, using- q = q, the spatial average of pruie is 
p/q = /, so that the shifted nu fields have zero spatial 
average. This shift is allowed and canonical because p is 
an integer. After this shift, the Hamiltonian becomes, up 
to a constant, 

'Hq = 

-t cos(A„(^i^ - + Ai-^) 



1 



+u nu -jjL mu- 



with jl = —{jls — p)/2. The Gauss' law constraint is 

changed to 



(F3) 



At this point, we can recover the analysis of Sec. Ill B 2 
by assuming the bosons are gapped, hu « 0, and treating 
the pure compact U{1)'^^^ gauge theory with the static 
source charges eu. Here, we will connect to the dual 
g-vortex theory directly. This is accomplished by per- 
forming a separate duality transformation for each of the 
g fractional boson fields. This is conveniently done, as 
elsewhere in the paper, in a path integral formulation. 
The appropriate world-line form of the partition func- 
tion, analogous to Eq. (2.5), but with the gauge fields 
included is 

Z= Y: fdA'l[S{A.jf)l[5(YBi]e-^^,iF4) 
where 



E- 



Jf -fl ( mj£-i I r 



2e2 

i{jf - Jt' + eur) ■ Al 



.B- 



s2 



+ E^^-(^x-^')« + tI^«i 



at 



(F5) 
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Here and in the remainder of this section, we wiU suppress 
3- vector indices and spatial coordinate labels when possi- 
ble, and use the usual vector notations (•, x) for 3- vector 
operations. Physically, = {nf ,Jg) is an integer-valued 
vector field representing the 3-current of fractional boson 
£, = {Aq,A^) is the 27r-periodic 3- vector potential for 
the i?*'* gauge symmetry. The is a 27rx integer- valued 
field defined on the links of the dual lattice, which, when 
summed over, gives a Villain potential for the f*'' gauge- 
field kinetic energy. Finally, f = (1,0,0) is the unit vec- 
tor in the time direction. 

To proceed, we solve the divergence constraints on the 
boson currents by defining vector potentials A^^ on the 
dual lattice, 



J = 



(A X A% 
27r 



(F6) 



and for convenience we find some fixed non-fluctuating 

g 

27rx integer fields B on the dual lattice satisfying 



with 



<j2 



= y- 



87r2e2 



(A X A^)i - 2'Kjl ( mu-- I T 



Ai+' + B;-A,Xi\ 



(F13) 



Next we "soften" the 27rx integer constraints on xti 
converting the partition sum to an integral, and adding 
terms in the action to weight integer values more heav- 
ily (alternatively, this can be done exactly by Poisson 
resummation, but we prefer to keep things simple): 



Z = JdA'dxelldl^Xa}! 



(F14) 



with 



S^„ = Si 



E 



[Vqy COS Al^ A cos Xae] ■ (F15) 



(A X B )i 
2tt 



(F7) 



Since J2e = 0> the background fields can (and should!) 

(, 

be chosen to satisfy B = 0. Inserting these into the 
partition function, one has 



{Al^,Bt} an \ I I 



(F8) 



ith 



y— 



87r2e2 



(A X A^)i - 27r/x mu 



-—Ai-[Ax{A'-A'+'+B 
27r ' ^ ^ 



B' 



§2 



(F9) 



At this point one may integrate out the A^ fields, which 
gives a delta-function constraining 

AxB^ = Ax{A^- A^+i + B^). (FIO) 
The general solution is 



At this point it is advantageous to perform a few manip- 
ulations in a row to bring the partition fimction into a 
simple form. Rather than show the results of the individ- 
ual steps, we simply give the procedure, to be executed 
in order: 

1. Decompose the background fields into longitudinal 
and transverse parts according to B^ = Aij,r]i + Bj^, 
where A • = 0, and one chooses J2e % = 0- 

2. Shift Xai Xae - i^ai + '&a,e+i + Vae, and integrate 
the partition function over the 'dae so introduced at 
every dual lattice site. 



3. Next shift z?^ dg, 

4. Finally, shift Af, - 



A,,^t. 



At this point, all explicit xt dependence will be elimi- 
nated from the action, and the integral over xi can be 
done to give an unimportant constant multiplier. The 
rescaled partition function is then 



-I 



(F16) 



with 



E 



(A X A% 



27r/i (rriie-- 
\ 1 



B,. 



AuX€, 



(Fll) 



where xe is a 27rx integer- valued scalar field. To preserve 
the constraint J2e^^ = 0> '^^ should impose J2iXe = 
(or at least a constant). The partition function becomes 



^= E HMEx- 

{Ai,xe} a \ e J 



(F12) 



■EtK 



at 



(F17) 



- E [^9" cos(A^i?£ - A^) zosipi - -dt+x - %)] • 



al 



We have nearly achieved the desired dual form. Indeed, 
due to the "Higgs" term, the only gapless gauge fluc- 
tuations are those with A^ = A equal for all £. The rel- 
ative fluctuations of the different A^ fields will be small 
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and can be neglected. There is in general, however, a 
mean value (average over configurations) for each of the 
fields, or more importantly the fluxes A x . This 
mean value is oscillatory and has zero spatial average, 
since both the ran — ^/l term has zero spatial average 
as docs A X = 27rei^f. Once can show by a consider- 
ation of the vortex band structure that, generically, the 
uniform vortex modes are essentially unaffected by this 
spatially oscillating flux with zero average. A careful in- 
spection shows, however, that the staggered potential Hs 
is crucial here. For the non-generic case Hs = 0, so that 
fi = p/2, one can show that the average flux is 

((A X A%) = 2np (^mu - , for /x, = 0. (F18) 

Although this has zero spatial average, this zero value is 



deceiving, since the background flux A x A^ differs on 
different plaquettcs by a multiple of 27r. Since the vor- 
tex kinetic energy is 27r-periodic, such a "zero average" 
flux has the effect of a uniform one, and drastically mod- 
ifles the vortex spectrum. In the generic situation with 
/is 0, however, the differences in fluxes between pla- 
quettes is truly non-trivial, and it is correct to neglect 
this oscillating average flux. 



In this case, we may take the continuum limit by set- 
ting A^ « A, e*"** — > |^|~"'^C^, adding potential terms that 
keep the magnitude of Q approximately constant, and 
assuming slow spatial variations of the Q fields. One ob- 
tains, after coarse-graining, an effective potential of the 
same form as Eq. (3.33) in the main text. 
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In fact, any of these operators can be composed with a (7(1) 
gauge transformation without any physical consequences. 
If we wish the gauge field to transform normally under 
the space-group operations of the PSG, this freedom re- 
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